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Abstract 

A measurement of the mass of the W boson is presented which is based on 
a sample of 5982 W — > ev decays observed in pp collisions at y/s = 1.8 TeV 
with the D0 detector during the 1992-1993 run. From a fit to the transverse 
mass spectrum, combined with measurements of the Z boson mass, the W 
boson mass is measured to be Mw = 80.350 ± 0.140 (stat.) ± 0.165 (syst.) ± 
0.160 (scale) GeV/c 2 . Detailed discussions of the determination of the abso- 
lute energy scale, the measured efficiencies, and all systematic uncertainties 
are presented. 
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I. INTRODUCTION 



Among electroweak measurables, the mass of the W boson Mw is of crucial importance. 
Along with the determination of the mass of the top quark |l]||] and in conjunction with 
other precisely determined quantities, including the mass of the Z boson M z the electroweak 
Standard Model [|J is constrained. This paper discusses the details of the first measurement 
of M\y by the D0 collaboration using data from the 1992-1993 running of the Fermilab 
Tevatron Collider. It includes essential calibrations which will be used in future D0 mea- 
surements. A first report of this measurement was published in Ref. 0]. 

An early success of the CERN pp collider was the discovery and measurement of the 
masses of both the W and Z bosons Table | gives a history of the published values of 

the direct measurements of Mw- The approach taken in this analysis is similar to those of 



the UA2 H and CDF [0,0 experiments. 

For the physics of W bosons, the electroweak measurables of interest are M w and sin 2 6w, 
where 9w is the weak mixing angle. Both can be measured precisely and can be predicted 



from the lowest order relations of the model 17 



M W = M Z cos 6 W (1) 
g 2 sm 2 6 w 

"em = : {*) 

47T 

V2 8M 2 / {6) 
Here, a EM is the fine structure constant, g is the gauge coupling associated with the SU{2) L 
gauge group, and is the Fermi coupling constant. The weak coupling, the electric charge, 

and the weak mixing angle are related by tan#iy = g'/g, where g' is the coupling of the 
U(l) gauge group. 

The standard set of measurable input parameters is the following: 

a EM = 1/(137.0359895 ± 0.0000061) (4) 

= 1.16639 (±0.00002) x 10~ 5 GeV~ 2 (5) 

M z = 91.1884 ± 0.0022 GeV/c 2 . (6) 



The fine structure constant is measured from the quantum Hall effect [18[; the Fermi coupling 



constant is measured from the muon lifetime 18], and Mz is measured directly by the 



combined LEP experiments |19] 



In order to confront the model beyond the lowest order, a self-consistent theoretical 
scheme for dealing with the effects of higher orders of perturbation theory is required. Dif- 
ferent theoretical prescriptions motivate particular definitions of the weak mixing angle. 
The determination of M\y and the ratio of neutral to charged current cross sections in deep 
inelastic neutrino scattering are most naturally interpreted in terms of the "on shell" scheme 
20fl in which the weak mixing angle is defined by Eq. ([[]) with both masses as measurable 



quantities. 

Within a given scheme, radiative corrections are most easily included through the use of 
a single measurable parameter, Ar, which is analogous to (g — 2) in quantum electrodynamic 
radiative corrections. Like (g—2), Ar can be determined experimentally and its measurement 
can be directly compared with its theoretical prediction. In the leading log approximation, 
it can be written in terms of Mw as 

A _ 1 7roW v2 , s 

G,Mi(M w /M z y\i - (M w /M z y\ ■ 1 ' 

Roughly 90% of the value of Ar is due to light quark loop corrections to a EM , while the 
balance is due to the embedded physics of heavy quarks and the Higgs boson. Any physics 
beyond the Higgs boson and known heavy quarks would also contribute to Ar. Prior to 
the measurement described here, the world average of M w (M w = 80.33 ± 0.170 GeV/c 2 ), 
M z HI, and Eq. (0), results in Ar = 0.0384 ± 0.0100, which is 3.8a from the tree level 



prediction. 

Because Ar is dominated by QED corrections, it is interesting to separate out those 
"residual" effects which are distinguishable from electrodynamic effects alone. Such a sep- 
aration isolates possible new physics as well as physics directly associated with the top 
quark and the Higgs boson. A prescription for doing this has been suggested by defining 
(Ar) rcs §nj, 



With a determination of Ar plus a separate evaluation of a EM (M|) p2| , (Ar) res can be 
extracted. Evaluation of this quantity makes a particularly economical probe of the Standard 
Model possible. 



A. Plan of the D0 Measurement 

The D0 detector is a calorimetric detector with nearly full kinematical coverage for elec- 
trons, hadrons, and muons. The inner tracking region does not include a magnetic field. 
Calibration of the electromagnetic, and by extension the hadronic, calorimeter was accom- 
plished by exposing calorimeter test modules to charged particle beams of known energies 
and compositions, as well as in situ decays of known particles. The D0 determination of 
M w relies on the determination of the mass ratio M w /M z and the subsequent scaling of 
this ratio by the precisely determined Mz from LEP |19| . This approach is similar to that 
of the UA2 experiment The significant advantage of determining Mw/Mz is that a 



number of systematic uncertainties cancel in the ratio. This paper addresses those in detail. 

The production and decay characteristics of W bosons in a pp collider present a variety 
of challenges which drive the analysis strategy. Because the statistical power of this mea- 
surement is at the level of less than 150 MeV/c 2 (< 0.2%), it is necessary to understand 
both the experimental and theoretical systematic uncertainties to a precision comparable to 
this level. The Z boson data are used for studying many of the experimental uncertainties, 
so the total uncertainty of an Mw determination is strongly coupled with the size of the Z 
boson data set. Hence, future determinations will gain in statistical and systematic precision 
with the sizes of both the W and Z boson data sets. 

Uncertainties in modeling the production of W and Z bosons present a different set 
of challenges which do not necessarily scale with the number of events. For example, at 
Tevatron energies roughly 80% of the annihilations which produce W bosons involve sea 
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quarks. Additionally, the substantial probability of gluon radiation from the initial state 
quarks results in significant transverse momentum of the W boson. Both of these theoretical 
issues involve uncertainties which complicate the simulation. 

For the Z boson the observables come from the reaction chain: 



where Hz is the hadronic recoil against the transverse motion of the Z boson. Both electrons 
are fully measured and the dielectron mass is determined from 



With the D0 detector, the electron decay mode of the W boson leads to the most precise 
mass determination. This is due to the cleaner signal and better resolution, as compared to 
the muon or tau decay modes. In this experiment, the single relevant channel is W — > ev e . 
The electrons are emitted with transverse momenta p e T of order 40 GeV/c and the neutrino 
is emitted with a comparable momentum, escaping without detection. This leaves a large 
component of missing energy $ in the event of which only the transverse component, $ T 
is determined. Therefore, the defining characteristics of W bosons are a high p T electron 
accompanied by significant $ T . 

The observable quantities for this measurement come from the W boson production and 
decay chain: 



where Hw is the hadronic recoil against the transverse motion of the W boson. Since a 
complete characterization of the neutrino 4-momentum is impossible, the only quantities 
directly measured are the electron momentum and the transverse momentum of the recoil 

Using these measurables, the two body kinematics of this decay provide at least two 
methods for measuring Mw- The transverse energy spectrum of the electron will exhibit a 



pp — > Z(— > e + + e ) + Hz(^ hadrons), 



(9) 




(10) 



pp — > W(— > e + v) + Hw{-^ hadrons), 



(11) 



9 



kinematical edge (the "Jacobian edge") at Mjy/2 for W bosons with transverse momentum 
p^ equal to zero. However, resolution effects and nonzero values of pYf smear the p e T and 
ifjrp spectra and therefore affect the use of the sharp edge as a measure of Mw- 

To control the systematic effects while retaining the highest statistical precision, the 
"transverse mass" is used to determine M w . It is defined as j23 



Ml = 2E e T E v T - Ip^p^ 

= 2£^(l-cos0 e „), (12) 
where p^ is the transverse momentum of the neutrino and 4> eu is the azimuthal angle between 



the electron and neutrino [2^]. The transverse mass also exhibits a Jacobian edge, but at 
the value of Mw and with much less sensitivity to p^ . Hence, precise determination of the 
location of this edge determines Mw- The effect of both the finite width T(jy) and p^ does 
distort the shape of the M? spectrum |25j . While the transverse momentum of the W boson 
is relatively low, peaking at approximately p^ ~ 5 GeV/c in this analysis, even this small 
amount can be significant. 



Equation [L2] shows that the necessary ingredients for determining Mt are p^, p£, and 
the angle between them. Among these, only p^ is determined directly. Since momentum 
transverse to the beam is conserved, a measured imbalance can be attributed to the neutrino. 
Therefore, in the absence of detector effects, the neutrino transverse energy is equal to the 
missing transverse energy and calculated from the measured p|? c and p%, 

p^ = $ T =-p™ c -Ft. (13) 

However, the reaction given in Eq. (|TTD does not fully describe the situation since energy 
measurement in a calorimeter includes other effects. Energy lost in detector cracks and 
inefficiencies can introduce biases in the magnitude and direction of the total energy. The 
interactions of the remaining spectator quarks in the proton and antiproton will add energy, 
as will noise and "pileup" due to the residual energy from multiple interactions. Designating 
these additional, non-recoil, luminosity dependent contributions the "underlying event", 
Ut{£>), the measured neutrino transverse momentum is given by 

10 



IPF + ViTiQ] -[*t(Q+ ft] ■ 



(14) 



The two quantities within the brackets are not distinguished from one another in the mea- 
surement but must be dealt with in the analysis. Figure [l] shows the kinematics of the W 
boson events. 

Since there is no analytic description of the transverse mass distribution, determination 
of Myy relied on modeling the transverse mass spectrum through a Monte Carlo simulation. 
The W boson mass was extracted by comparing the measured distribution in transverse mass 
to the Monte Carlo distribution generated for different W boson masses. The simulation 
relied on experimental data as much as possible and used Z boson events, not only to set the 
energy scale, but also to understand the electron energy resolution, the energy underlying 
the W boson, and the scale in p^. . In the simulation W bosons were generated with a 
relativistic Breit-Wigner line shape that was skewed by the mass dependence of the parton 
luminosities. The longitudinal and transverse momentum spectrum were given by a double 
differential distribution calculated to next-to-leading order. 

The decay products and the W boson recoil system were traced through the simulated 
detector with resolution smearing. Minimum bias events (collisions which are recorded with 
little or no trigger bias) mimic the debris in the event produced by the spectator quarks 
and pileup associated with multiple interactions. The minimum bias events also properly 
included residual energy from previous crossings. The generated spectra in transverse mass 
for different values of My/ were compared to the measured spectra by a maximum likelihood 
method, and the best fit value of the mass obtained. 

The measurements reported in this paper are: Myy, as determined from fits to the 
distribution and fits to the p\. and p? distributions and Mw/Mz- In addition Ar and (Ar) res 
are determined. 

The paper is organized as follows: 

• Section II: a brief description of the detector; 
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• Section III: data collection, reconstruction, the corrections applied to the data and the 
selection of the final sample; 

• Section IV: determination of the parameters used in the Monte Carlo simulation; 

• Section V: Monte Carlo simulation; 

• Section VI: results of the fits; 

• Section VII: effects of systematic errors due to the parameters determined in Section 
IV and the assumptions described in Section V; 

• Section VIII: consistency checks; 

• Section IX: conclusions; 

• Appendix A: W Boson and Z Boson Production Model; 

• Appendix B: Bremsstrahlung; and 

• Appendix C: Mean Number of Interactions. 

II. D0 DETECTOR AND TRIGGER SYSTEM 

The data collected for this measurement were taken during the exposure of the D0 
detector to collisions of protons and antiprotons at a center of mass energy of ^/s = 1800 GeV 
in the 1992-1993 running period of the Fermilab Tevatron collider. This was the first beam 
exposure of this experiment and the total luminosity accumulated was 12.8 pb _1 . The 
average instantaneous luminosity was C = 3.4 x 10 30 cm~ 2 s~ 1 , which corresponded to less 
than 1 collision per crossing which occurs for £~ 6 x 10 30 cm~ 2 s _1 

The D0 detector was designed to study a variety of high transverse momentum physics 
topics and has been described in detail elsewhere |2l| . The detector has nearly full accep- 
tance for electrons, photons, and muons and measures jets, electromagnetic (EM) showers, 
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and ifjrp with good resolution f27| . The detector consists of three major subsystems: a 
tracking system, uranium-liquid argon calorimeters, and a muon toroidal spectrometer. The 
components of the detector which are most relevant to this analysis are briefly described 
below. 



A. Tracking System 

The tracking system was used to reconstruct charged particle tracks over the region 



I77 1 < 3.2 and to reconstruct the interaction vertex of the event. It consists of four 
subsystems: a drift chamber surrounding the vertex region (VTX), a transition radiation 
detector (TRD), a central drift chamber (CDC) and two forward drift chambers (FDC). The 
VTX, TRD, and CDC cover the large angle region and are oriented parallel to the beam axis. 
The FDC's cover the small angle region and are oriented perpendicular to the beam axis. In 
addition to the r<p measurement of hits in the CDC, delay lines were used to measure track 
hit locations in the z direction. The TRD provides an independent identification of electrons, 
in addition to that provided by the calorimeters, allowing enhanced hadron rejection. 



B. Calorimeters 

The calorimeter system consists of one central (CC) and two end (EC) calorimeters which 
measure the energy flow in the event over a pseudorapidity range \rj\ < 4.2. The calorimeters 
are enclosed in three separate cryostats which surround the tracking system. They each have 
an electromagnetic, a fine hadronic (FH), and a coarse hadronic (CH) section. Liquid argon 
is employed as the active medium and uranium is the absorber material in the EM and 
FH sections and copper (steel) is the absorber in the CH section for the CC (EC). The 
inter-cryostat region (ICR) is instrumented with scintillator tile detectors which are located 
in the space between the EC and CC cryostats. These detectors were used to improve the 
energy measurement of jets that straddle two calorimeters. 
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The calorimeters are arranged in a cylindrical geometry with each EM section being 
divided into four longitudinal readout layers, for a total depth of 21 radiation lengths. A 
projective tower arrangement for readout points toward the interaction region. The hadronic 
sections are 7-9 nuclear interaction lengths deep and are divided into four (CC) or five (EC) 
longitudinal readout layers. The transverse segmentation of the calorimeters is 0.1 x 0.1 in 
A77 x A0, except in the third layer of the EM calorimeter which is at shower maximum, 
where it is 0.05 x 0.05 in Arj x A0. Measured resolutions will be discussed below. 

C. Muon Spectrometer 

The muon spectrometer provides identification and momentum determination for muons. 
It surrounds the calorimeters and consists of planes of proportional drift tubes which sur- 
round magnetized iron toroids and covers a region |?7| < 3.3. There is one layer of propor- 
tional tubes on the inner face of the magnet and two layers, separated by « lm, outside the 
magnet. The material in the calorimeter and iron toroids is about 12 interaction lengths 
thick, making hadronic punch-through to the outer two layers negligible. The muon mo- 
mentum resolution is a(l/p) = 0.18(p — 2)/p 2 © 0.008 (momentum p in GeV/c). 

D. Triggers 

The Tevatron beam crossings occurred every 3.5 /is. For a pp total cross section at 
a/s = 1.8 TeV of approximately 70 mb |f29fl , there is an interaction rate of ~ 200 kHz at a 
typical instantaneous luminosity of 3 • 10 30 cm~ 2 s _1 . In order to record events at ~ 2 Hz, 
three stages of hardware and software triggers were used. 

To indicate the presence of a collision within the detector and to calculate a fast approx- 
imation to the vertex position, radial scintillation hodoscope arrays are positioned in the 
forward directions subtending angles of 2.3 < \r]\ < 3.9. To pass the level (L0) hardware 
trigger, coincident hits in these counters were required on both sides of the interaction region, 
signaling an inelastic collision within the detector volume and also providing an estimate of 
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the z position of the interaction vertex. This trigger provided the minimum bias data set 
used in this analysis. The trigger rate depended on the luminosity and for the data analyzed 
here, was typically 90 kHz. 

Events passing the L0 trigger were then passed to the level 1 (LI) hardware trigger. Here 
a decision was made based on the fast analog sums of all the EM layer calorimeter signals 
which represent the energies in trigger towers. These towers were segmented as 0.2 x 0.2 in 
Ar] x A0 with coverage extending to \rj\ < 3.2. The LI electronics restricted the maximum 
trigger rate to ~ 200 Hz and decisions made by it and by the L0 trigger were made between 
beam crossings. 

The final stage of triggering was the level 2 (L2) software trigger which ran on a farm 
of 48 VAXstation 4000 M60 processors. The typical processing time for an event in L2 was 
350 /is, resulting in an average deadtime of ~ 2%. The full segmentation of the calorimeter 
was available at this trigger level and a full event reconstruction was done, albeit with 
simplified algorithms and coarser segmentation. There were 32 different LI components and 
128 different L2 components which could be constructed and prescaled at different rates. 
Each logical combination targeted specific physics for given accelerator conditions. These 
data sets were written out to corresponding output data streams. 

1. Trigger requirements for the W boson and Z boson data sets 

For the determination of the W boson mass, the electron decay modes were required 
for the selection of both W and Z bosons. While the characteristics of both are similar— 
the presence of high transverse momentum electrons — the different rates and backgrounds 
require distinct selection criteria. Because of the presence of a neutrino in W boson decay 
events, a minimal $ T requirement was used in the selection of W boson candidates. 

To select electrons, the LI trigger required the transverse energy in the EM layers of a 
trigger tower to be above a preselected threshold. For the selection of W boson events, the 
LI trigger required at least one EM trigger tower above 10 GeV. For Z boson events, at 



15 



least two EM trigger towers with Et > 7 GeV were required. 

The L2 electron algorithm used the full segmentation of the EM calorimeter to measure 



the energy deposited by the EM shower and is described in detail in Ref. [[Z7| . Trigger towers 
above threshold were used as seeds to form energy clusters which included all calorimeter 
cells in the four EM layers and the first FH layer in a window Arj x A0 = 0.3 x 0.3. For 
the selection of W events, an energy cluster with Et > 20 GeV was required by the L2 
filter. For Z boson event candidates, two energy clusters, each with Et > 10 GeV, were 
required by the L2 filter. Transverse and longitudinal shape requirements as well as isolation 
requirements were also placed on the energy cluster for the selection of W boson events. 

The $ T in the event was calculated in the L2 trigger and was required to be above 20 
GeV for the W boson event selection. It was computed using the vector sum of the Et of 
all the cells in the calorimeter and the ICD with respect to the z position of the interaction 
vertex, determined by L0. Prescaled triggers of W boson events were also recorded without 
the Ifirp trigger requirement in order to study efficiencies and biases. Table [TT] lists a summary 
of the trigger requirements used in the selection of the W boson and Z boson data samples. 



2. Main Ring veto 

The 150 GeV/c conventional accelerator (main ring) passed through the coarse hadronic 
part of the D0 calorimeters. Halo particles accompanying the circulating beam can deposit 
energy in the calorimeter and corrupt measurements both at the trigger and offline recon- 
struction levels. Such unrelated energy depositions in a localized part of the detector will 
affect the $ T determination and therefore considerable care was required in the utilization 
of triggers taken while the Main Ring beam passed through the detector. A veto gated on 
the injection period of the main ring cycle (the first 0.4s of the 2.4s cycle) was used in some 
of the L2 filters to avoid any adverse effects. For the mass analysis, all Z — > ee events were 
used; Z — > ee candidates recorded during the veto of the main-ring period were excluded 
from resolution studies. No W boson events were taken from triggers occuring during the 
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main ring cycle. 



III. W AND Z BOSON DATA SAMPLE: EVENT SELECTION AND DATA 

PROCESSING 

A. Offline Data Processing and Candidate Event Selection 

Once data were written to tape, the digitized information was converted using an initial 
calibration. The initial calibration of the CDC was based on measurements using cosmic 
ray data. The calibration of the delay lines of the CDC, measuring the z coordinate of the 
hits, was derived from pulser measurements on the bench combined with cosmic ray data. 



The initial calibration of the calorimeter was obtained from test beam measurements |26 
Corrections to the calorimeter calibration were required, including corrections for an adjust- 
ment to the operating voltage of the calorimeter, corrections to the sampling weights, and 
to the gains of individual calorimeter cells. In addition, a correction due to a difference in 



liquid argon temperature at the test beam and D0 was applied [jJlJ. It should be noted 
that none of the central EM calorimeter modules that were tested were installed in the final 
calorimeter. 

The azimuthal uniformity of the central electromagnetic calorimeter was determined 
using approximately 3.5 million triggers from an inclusive electron data sample [JT|. By 



equalizing the event rate above a 13 GeV threshold for each calorimeter module, relative 
calibration constants were determined to an accuracy of 0.5%, assuming that the observed 
variations were instrumental in origin. These relative calibration constants showed a 
variation in the response between different modules with a maximum difference of 5%. The 
variations were dependent on which of the 32 EM modules was struck by the electron, and 
not by other features of the calorimeter such as a variation in the amount of material in the 
tracking detector. All of the above corrections to the energy are propagated into the $ T 
calculation. 
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1. EM Clustering 



Electrons and photons were reconstructed as energy clusters in the EM and first FH 
section of the calorimeter. Towers were defined by adding the energy measured by the 
calorimeter in all four EM layers plus the first FH layer for cells within 0.1 x 0.1 in A 77 x A0. 
Towers were grouped together with their adjacent neighbors, provided their energies are 
above 50 MeV. Clusters of adjacent calorimeter towers with significant energy depositions 
were then formed using a nearest neighbor clustering algorithm |H2J . 

The observed energy of the EM cluster is given by 

E meaa = 5 TB + Cj2PiSi (15) 

i 

where 5 TB is an offset in the energy response due to energy loss in the material upstream 
of the test calorimeter. C is the conversion constant from the digital signal to energy and 
(3i are the sampling weights for the i th layer with energy deposition (in ADC counts) 
The sum runs over all five layers in the EM calorimeter which contribute to the EM cluster. 
Both fa and 5 TB were determined from test beam measurements using electron beams over 
a broad range of energy and rapidity. From these test beam measurements, the offset was 
determined to be 5 TB = 347 MeV. This level of energy scale determination, based on the 
test beam measurements and taking into account corrections due to the transfer of the 
calibration from the test beam to D0, resulted in an energy scale approximately 5% lower 
than the nominal energy scale for the central calorimeter, as observed from the measured Z 
boson mass. Both the final overall energy scale and the offset were re-computed in situ, as 
will be discussed in the next section. 



2. Electron Identification 
Electrons were identified by a combination of topological and kinematic identifiers as 



described in [^7]]. The main electron identification requirements are below . 
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The cluster shower shape can be characterized by both longitudinal and lateral energy 
depositions. The fraction of the cluster energy which is deposited in the EM calorimeter is 
defined as Jem- Since charged hadrons deposit less than « 10% of their energy in the EM 
calorimeter, /em provides a powerful discriminant. Jem was required to be greater than 



0. 



The electron candidate cluster was required to have a topology, both longitudinal and 
lateral, which was consistent with that of electrons from a detailed GEANT Monte Carlo 
simulation |33] which was extensively compared to test beam measurements. A covariance 



matrix of 41 observables was defined to characterize an electron shower [27j. A x 2 parameter 
was defined to measure the consistency of the shower with that expected for electrons. For 
central electrons a x 2 < 100 was required; for electrons in the end calorimeters x 2 < 200 was 
imposed. These requirements were ~ 94% efficient. A rejection factor of about 4 against 
EM clusters that are not due to single electrons was achieved. 

Electrons from W and Z boson decays tend to be isolated from other particles in the 
electromagnetic calorimeter. An offline isolation requirement was used which is defined as 

, E to t{Ri) — Eem\Rv 

where E tot is the total energy in cone of radius R4 = 0.4 in r\ x 0-space. E E m is the EM 
energy in a cone of radius R2 = 0.2. An isolation requirement of /; so < 0.15 was placed on 
the cluster energies for electrons from both W and Z boson event candidates. 

An important source of background for electrons is photons from tt° or 77 meson decays 
which are adjacent to unrelated tracks. This background was reduced by requiring that a 
track from a charged particle in the tracking detector be consistent with the position of the 
cluster in the calorimeter. To qualify as a match, the shower centroid was required to link 
with a reconstructed track with significance <Jtrk, 



®trk = 



+ (17) 



\ \cr RAip ) 

where RAtp and Az are the spatial mismatches between the track projection and cluster 
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position in the and z directions, respectively, and ora^ and a z are the associated ex- 
perimental resolutions. For the data set used in this analysis a t rk < 10 was imposed for 
CC electron candidates. This cut was used to minimize bias and results in an efficiency of 
> 98% pl. 



3. Electron Direction Determination 

The optimum resolution in the electron polar angle is obtained using the z position of 
the electron cluster as obtained from the calorimeter information and the z position of the 
center-of-gravity (cog) of the CDC track. These two points thus define the polar angle 
of the electron. The position of an electron in the calorimeter was determined from the 
energy depositions in the third EM layer of the shower using a log(E) energy weighting 
algorithm ||35|| . The parameters of the algorithm were determined using both test beam and 
Monte Carlo data. Further study with collider data demonstrated the need to remove a 
residual bias in the z position of the cog of the track. This was accomplished in situ using 
an inclusive muon data set which demonstrated that the z position of the cog of the CDC 
track was biased. 

The true and measured z position of the cog of the track are related linearly by 

Zfrue = OCDC 2 meas + PcDC- (18) 

With the muon data sample, the scale factor ckcdc was determined by defining a track using 
the cog of the muon track, as reconstructed in the first layer of muon spectrometer before 
the toroidal magnet, and the vertex z position. By comparing the expected and measured 
z positions the scale factor was determined to be «cdc = 0.988 ± 0.002, where the error 
is the combination of a small statistical component (±0.0003) and the following systematic 
components: the observed variations in scale factor for different azimuthal regions of the 
detector (0.001), observed variations for different polar angles of the muon tracks (0.001), 
muon chamber alignment (0.0003), and different methods to extract «cdc (0.0004). The 
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offset /^cdc is consistent with zero. In this analysis cxcdc = 0.988 and (3cdc = were 
used. The scale factor has been confirmed using cosmic ray muon data. The trajectory of 
cosmic ray muons traversing the full detector was reconstructed using the non-magnetic inner 
volume of the spectrometer. As before, the expected CDC track positions were compared to 
their measured positions and the scaling behavior of the cog of the CDC track was, within 
its uncertainty, confirmed. 

To verify the consistency of the electron angle determination, Z — > ee events were stud- 
ied. Given the electron cluster position, their track intersections with the beamline were 
reconstructed from the Z — > ee decay. These intersections in general do not coincide. After 
applying the correction to the CDC cog, the width of the distribution in the difference in z 
positions of the two intersections was tracked by varying the calorimeter electron z position. 
The resultant width of the distribution was at a minimum without applying any correction 
to the electron calorimeter position, showing the internal consistency of the procedure. 

4- Measurement of $ T 

The total missing transverse energy in the event is calculated by summing over all 
calorimeter and ICD cells 

$ T =-Y / E l sm6 i u i = -Y / E l T . (19) 

i i 

Here Ui is a unit vector from the event vertex to the cog of the i-th calorimeter cell, Ei the 
energy in the i-th calorimeter cell, and 9i is the polar angle given by the event vertex and 
the cog of the i-th calorimeter cell. 

The nominal event vertex was determined using all tracks in the CDC. The transverse 
momentum of the electron was calculated using the total energy and direction of the cluster 
in the calorimeter. Since, the electron direction, as computed above, may not intersect the 
nominal event vertex, a recalculation of the $ T was done by using the vertex obtained from 
the electron alone. For Z — > ee events, the event vertex and the electron polar angles were 
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determined using a constrained fit of the measured variables of the two electron directions. 
The missing transverse energy used in this analysis was based on calorimetric information 
alone and was not corrected for possible muons in the event. 

B. Final W — > ev and Z — ► ee Data Sample 

After electron identification and calculation of the missing transverse energy, the final 
W boson and Z boson candidate samples were subjected to the following selection criteria: 

Fiducial requirements were placed on the electrons in the W boson candidate events to 
select central electrons: < 12. Here, is an integer index for the calorimeter tower 
containing the most energetic cell of the electron cluster in the third EM layer. It is equal to 
10 x 77 for particles which originate at z — 0. In order to ensure no energy leakage into the 
uninstrumented regions within modules, electrons were restricted from the readout edges in 
(f> by requiring that their impact to be within the central 80% of each module. An additional 
requirement was imposed that no event have a jet in which the fraction of energy in the 
CH section of the calorimeter exceeds 0.4. This eliminated events with spurious energy 
depositions from the Main Ring. 

The kinematic and fiducial requirements that defined the W boson candidate sample are: 

• K\ < 12 

• Ej, > 25 GeV 

• $ T > 25 GeV 

• Vt < 30 GeV/c 

This resulted in a sample of 7262 events. Additionally, including a transverse mass cut of 
Mt < 110 GeV/c 2 left 7234 candidates. 

Z boson candidate events were accepted with the requirements: 

• \i*\ < 12 or \i e v \ > 15 
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• E^ 2 > 25 GeV . 

As in the W boson sample, the module boundary edge cut was made for CC electrons only. 
For the final Z data sample, only events with both electrons in the CC (|^| < 12) were 
used. This resulted in a sample of 395 candidate events with both electrons in the CC. A Z 
boson mass cut which eliminated events outside a window of 70 < m ee < 110 GeV/c 2 left 
366 candidates. For some studies, events with one electron in the forward region (\i^\ > 15) 
were included. For resolution studies, Z — > ee events had the additional requirement that 
events were not accepted when taken during the Main Ring cycle. 

Table [TTT] lists the final number of events in the samples. Figure 0(a) shows the transverse 
mass distribution of the central W boson candidate events, before the transverse mass cut 
and 0(b) shows the invariant mass distribution of the central Z boson candidate events. 
Neither distribution is corrected for the electron energy scale determined in situ (see Section 
IV). 

C. Data Samples Used in the Analysis 

There were five primary data samples which are used in this analysis. 

W — > ev sample: A sample of 7262 W — > ev candidates (prior to fitting and transverse 
mass cut) provided the main data sample used to measure the W boson mass . 

Z — > e + e~ sample: A sample of 395 central Z —>■ e + e~ candidates (prior to the Z bo- 
son mass cut) was used along with the W — > ev sample to measure the W boson 
mass. A slightly enlarged sample was used in the determination of detector response 
parameters. 

Minimum Bias sample: A sample of approximately 50,000 triggers, taken at various lumi- 
nosities, was used for modeling the underlying event. 

J/i/; —> e + e~ sample: A data set of approximately 50 observed J/ip candidates was used 
in the electron energy scale determination. 
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7T° — > 77 — > e + e e + e sample: A data set of approximately 2500 observed 7r° candidates 
was used in the electron energy scale determination. 



IV. DETERMINATION OF PARAMETERS IN MONTE CARLO SIMULATION 

The extraction of Mw relied on an accurate and fast Monte Carlo simulation. The 
details of the physics model used in the simulation will be discussed in the next section. 
However, many parameters such as calorimeter response, efficiencies, and resolutions, were 
input to the simulation and were derived from the data. The focus of this section is a 
detailed description of how these parameters were determined. The use of these parameters 
in the simulation appears in Section V and the systematic errors on Mw due to uncertainties 
inherent in these parameters is described in Section VII. 



A. Electron Energy Scale 

All calorimetric measurements rely on the determination of the overall energy scale using 
particles of known momentum and/or on the reconstruction of the mass of well known 
particles. Both techniques have been used to calibrate the D0 calorimeter. Since the 
absolute energy scale of the EM calorimeter was not known to the required precision, the 
ratio of the measured W boson and Z boson masses and the world average Z boson mass 
were used to determine the W boson mass. A number of systematic effects, common to both 
measurements, cancel in the ratio. Most notably, as will be discussed in more detail below, 
the ratio was, to first order, insensitive to the absolute energy scale. 

The initial calibration of the calorimeter was provided by transferring the calibration from 
a test beam to the collider detector, as discussed in Section [III A| [[36|| . An important result 



of these test beam measurements was the demonstration that the EM calorimeter is linear 
to better than 0.5% for electron energies exceeding 10 GeV. To complete the establishment 
of the energy scale with the desired precision, it was necessary to determine to what extent 
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a possible offset in the energy response, as opposed to only a scale factor, was responsible 
for the deviation of the ratio M% /M^ EP from unity. 

A strategy for establishing the final energy scale and possible offset in the response was 
implemented. Inherent to this program is the assumption that the measured energy E meas 
is related to the true energy, E tTue , by a scale a and offset 5: 

E meas = aE tTac + 5. (20) 

Then, for a two body decay when 5 <C (Ei + E 2 ), the measured invariant mass of the decay 
products m mcas is related to the true mass m truc by 

m mcas « am tIUC + 5xf . (21) 

Here, / is a parameter that depends on the kinematics of the decay and is given by 

= (ff r as + jg mcas )(1 _ COS7) 

where £™| as are the measured energies of the two decay products and 7 the opening angle 
between them. When 8 is small, / is nearly equal to dm mea,s /dS. Hence, sensitivities to 8 
can be different, depending on /. 

Consequently, the dependence of the measured ratio of the W boson to Z boson masses 
on a, 8 can be estimated from the relation 

M w (a,8) 



M z (a,8) 
M w 



M 



z 



true 



1+ 6 f w M z -f z M w 



(23) 



a M z ■ M w 

Here, fw and fz correspond to average values of / for the W and Z bosons, respectively. 
Note that the determination of Mw from this ratio is insensitive to a if 5 — 0, and that the 
correction due to a non- vanishing value for 8 is strongly suppressed due to the fact that the 
W and Z boson masses are nearly equal. 

The values of a and 8 were determined from the analysis of collider events containing 
two-body decays for which 772 truc is known from other measurements. The three decays used 
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are the Z — > ee decays, measurements of 7r° — > 77 — > 4e states, and J/ip — > ee states. These 
three decays probe a useful range in /. The reference mass values used as benchmarks 
are: M z = 91.1884 ± 0.0022 GeV/c 2 0, m J/4) = 3.09688 ± 0.00004 GeV/c 2 0, and 
7?V) = 134.976 ± 0.0006 MeV/c 2 [|TJ. 

Z — > ee Analysis: The strongest constraint on the energy scale uncertainty comes from 
the Z boson data. The fact that electrons from Z boson decays are not monochromatic is 
exploited by studying the invariant mass distribution as a function of the variable fz- Small 
values of fz correspond to the decay of highly boosted Z bosons with, on average, higher 
energies. The dependence of the observed Z boson mass as function of fz, shown in Fig. || 
thus directly translates into a constraint on the energy scale and offset. This analysis was 
based on Z boson events with both electrons in the CC which were required to pass the 
same selection criteria as the final Z boson sample, except that Et > 10 GeV was required 
for both of the electrons. The data were binned in fz and the distribution in m™ eas was 
fit using a convolution of the Z boson Breit-Wigner resonance with a Gaussian resolution 
function. Using the standard Monte Carlo generator (described in the next section), sample 
distributions in m ee were generated under different assumptions for a and 5. A \ 2 comparison 
was performed between the data and the Monte Carlo and the lcr constraint on a and S 
from the Z boson data, shown as the solid line in Fig. was determined. 

7T° Analysis: ir° mesons were observed through their two photon decay and subsequent 
conversion to unresolved e + e~ pairs. There is a 10% probability for each photon to convert 
in front of the CDC, so that when both photons convert the dE/dx can be measured in the 
drift chamber and a strategy for the identification of 7r° decays is possible. The identification 
requirement was that one electromagnetic cluster be observed with two doubly ionizing tracks 
pointing to it. The diphoton opening angle, 8, was calculated from the center of gravity of 
those two tracks and the measured vertex of the event. In this way, an approximation of 
the mass was calculated as 

m^ = (£ clus )sin^, (24) 
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where E c \ us is the cluster energy which is equal to the sum of the photon energies since 
the photons are not resolved in the calorimeter. This strategy assumes a symmetric decay. 
Figure shows the signal and background, the latter determined from a single-conversion 
control sample. The invariant mass spectrum of the background-subtracted signal compares 
well with a Monte Carlo simulation shown as the solid line in Fig. ^. The measured mass is 
135.4 ± 10.0 MeV/c 2 . The sensitivity to the energy scale and offset is determined by 
varying both parameters in the Monte Carlo simulation and performing a x 2 fit to the data. 
Since the response given in Eq. (|2~0"D is the response per electron, the offset in response for 
the 7T° is 5 n o = 45. This procedure maps out an allowed region in the (a, <5)-plane shown as 
the dashed line in Fig. [|. 

J /if) Analysis: A sample of J /if) — *> ee events was also used in the EM energy scale 
determination. The data were collected in a set of special runs and had an integrated 
luminosity of ~ 100 nb _1 . The LI trigger required two EM triggers towers above a 2.5 
GeV threshold with less than 1 GeV in the corresponding hadronic towers. At L2, two EM 
energy clusters above 3 GeV were required and the isolation fraction was required to be 
/i SO < 0.4. Since the major background is due to n° — > 77 and r\ — > 77 decays in which one 
of the photons converts before it reaches the tracking chamber, the track associated with the 
electron cluster was required to have an energy deposition in the tracking chamber of less 
than 1.5 times the energy deposition of a minimum ionizing particle (MIP). In addition a 
cut was placed on the width of the cluster, defined as the average weighted distance of each 
cell of the cluster from its center. The weights are the same as those used in the position 
finding algorithm. The opening angle between the two electrons was determined from the 
event vertex and the cluster positions in the calorimeter. Figure |7| shows a clear J /if) — > ee 
signal above background. The background mass distribution for the J / if) signal was obtained 
independently by pairing EM energy clusters in the calorimeter in which at least one of the 
EM clusters had no associated track. If there was an associated track, it was required to 
have an energy deposition greater than 1.5 MIP. The remaining requirements imposed on 
the EM energy clusters were the same as in the analysis of the ee events. The mass value 
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fit for these data is rrij/^ = 3.032 ± 0.035 ± 0.190 GeV/c 2 , where the first error is statistical 
and the second is systematic. 



The model used for comparison to the mm distribution was an ISAJET |37|] based sim- 
ulation for bb production and subsequent decay to J/ip followed by a GEANT detector sim- 
ulation. Since there is additional jet activity close to the electrons in J/ip decays from b 
quarks, the contribution from the underlying event energy was evaluated using the simula- 
tion. Two classes of events were generated consisting of only the e + e~ pair of the J/ip decay 
and events corresponding to the full pp collision. The difference between the fully simulated 
and the reconstructed mass was 80 MeV/c 2 . This difference was applied as a corection and 
an uncertainty of 100% was assigned to this correction. Figure ^ shows the constraint on a, 
8 from the J/ip analysis indicated by the dotted line. 

Underlying Event Contribution: Background energy and noise can contribute to the 
measurements of electron energies. The different environments for Mz, m n o, and mj/^ final 
states led to different corrections for each. Monte Carlo studies specialized to the scale 
analysis plus the understanding of instrumental effects lead to corrections for Mz, m n o, and 
mj/4 of 0.17 ± 0.05, 0.30 ± 0.10, and 0.08 ± 0.08 GeV/c 2 , respectively. The uncertainties 
on these measurements form the dominant uncertainties in the determination of a and 8 for 
the 7T° and J / ip analyses. 

Combined Analysis: The data from the three samples are combined by adding the x 2 
distributions. For the combined x 2 the minimum value is x 2 — 53.8 for 58 degrees of freedom, 
with a best fit of 8 = -0.158±0.016 GeV and a = 0.9514±0.0018. This is consistent with the 
result obtained from the Z boson data only, but with substantially reduced errors. Figure |8| 
shows an enlargement of the region where the contours from the three data samples overlap. 
The shaded area is the contour obtained from the combined analysis for a unit change in 
the combined x 2 ■ 

The main contributions to the systematic uncertainties were the underlying event cor- 
rection and possible non-linearities in the energy response. Varying the underlying event 
by the errors quoted above changes the value of 8 for which the combined x 2 is minimized 
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by ±2 MeV when varying the underlying event contribution to the J /if) and by ±30 MeV 
when varying the contribution to the tt°. In addition the calorimeter 7r° response was varied 
taking S w o = 35, rather than AS as discussed above. This decreased the best value for 5 by 
52 MeV. The dominant uncertainty comes from a possible nonlinearity of the calorimeter 
and has been addressed by studying test beam data. The test beam data permitted a small 
nonlinear response of the EM calorimeter and was parametrized by including a quadratic 
term in the energy response of Eq. (|20|) , which was constrained by the test beam data to not 
exceed 1 part in 10 4 . Allowing for a nonlinear response characterized by such a quadratic 
term and repeating the above analyses results in an allowed region in a, 5 indicated by the 
dotted line in Fig. ||. The result is to decrease 5 by 200 MeV. The energy scale parameters 
and their uncertainties are thus 

5 = -0.158 ± 0.016±g;5i GeV, (25) 
a = 0.9514 ± 0.0018lg;^ (26) 

where the first error is statistical and the second systematic. The effect of a possible 
quadratic response term was included as the asymmetric contribution to the overall un- 
certainty shown on 5. 

The result described in this section constitutes the calibration of the central EM calorime- 
ter after the initial calibration based on a transfer of the calibration from the test beam given 
by Eq. fll"5|) . In practice, inserting the offset as defined in Eq. fli~5|) into Eq. (|20|) leads to 

E mcas = 5 TB + 5 + aE true (27) 

demonstrating that the in situ determination of 5 amounts to a redetermination of 5 TB . 
Combining Eq. (|27|) with Eq. (|l"5|) leads to 

£ truc = -- + - J £l3 i S i . (28) 
a a 

Using the redetermined values of a and 5, an overall offset of —5 /a = 158/0.9514 = 166 MeV 
was observed consistent with the average energy loss by electrons in the material before the 
calorimeter which was predicted by GEANT Monte Carlo simulation studies. 
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B. Hadronic Energy Scale 



The scale of the measured recoil momentum differs from the electron energy scale because 
the recoil measurement also includes energy from hadronic showers and suffers from the 
loss of energy in uninstrumented regions of the calorimeter. The response of the hadronic 
calorimeter relative to the response of the electromagnetic calorimeter was determined from 
Z —>■ ee events. In Z — > ee events the transverse momentum of the Z boson p^ can be 
obtained from either the measurement of the transverse momenta of the two electrons p^ e or 
from the recoil activity in the event ~Pt c - The latter was the way in which p^ was measured. 
To minimize the effects of the energy resolution in the determination of the hadronic energy 
scale relative to the electromagnetic energy scale, the momentum imbalance was measured 
with respect to the (rj, £ )-coordinate system |L2] . The rj axis is defined as the bisector of the 



two electron transverse directions. In the plane of the electrons, the axis orthogonal to the 
rj axis is the £ axis (see Fig. The rj imbalance is then defined as 

r]imb =Pt-V + Pt c ■ V (29) 

with fj a unit vector along the r\ axis. If the electromagnetic and hadronic responses are 
equal, ?7i m b is zero. Since the positive 77 axis is always in the direction of p^, any systematic 
bias in the measurement of pj? c will manifest itself as a bias in r/i m b- If the difference is due 
only to a scale, then the relationship between the two quantities can be characterized by a 
proportionality constant k. 

The determination of the hadronic energy scale factor requires selection of Z — > ee events 
with the same event topology as W — > ev events. Z — > ee events were selected with at least 
one electron in the central calorimeter. An additional cut to eliminate events which occur 
during the Main Ring cycle was imposed to ensure that no spurious calorimetric depositions 
affecting the measurement of the hadronic recoil were present. As a consistency check, Z 
boson events with both electrons in the central calorimeter have been used and a consistent 
result for the hadronic energy scale was obtained. Three related determinations of the 
hadronic energy response relative to the electromagnetic response have been carried out: 
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1. The primary method of obtaining the calorimeter response used was the measurement 
of the rj imbalance as function of \p^ e ■ fj\, as shown in Fig. |T0|(a). A least squares fit 
yields \pj? c ■ rj\ — k \p T e ■ fj\, with k = 0.83 ± 0.03. The offset in response, obtained 



from the intercept of the fit with the ordinate in Fig. |10|(a), was measured to be 
—0.17 ± 0.24 GeV. This result is consistent with zero. Figure 0(b) shows the i] 
imbalance for k = 0.83. The distribution is well described by a Gaussian distribution, 
centered at zero, with a width of 4.2 GeV. 



2. A second, very similar approach to fixing the scale k of the recoil system with respect 
to the dielectron system was the measurement of \p^ c ■ r)\ as function of \p T e ■ fj\, as 



shown in Fig. [11]. The linear dependence shows that, over the range of interest to 
the W boson mass measurement, the hadronic recoil is related to the electromagnetic 
energy by a simple scale factor. The scale n was determined by a least squares fit to 
the data, where the errors on p T ec ■ fj and p T e ■ fj have been determined using the known 
detector resolutions. This method gives k = 0.84 ± 0.03. The offset in response is 
0.06 ± 0.25 GeV, consistent with zero. It should be noted that the contribution from 
the underlying event ut does not affect the determination of k since it is distributed 
randomly with respect to the fj direction. 

3. The hadronic energy scale k was also determined using a third method which yielded 
both the hadronic energy scale and the magnitude of the underlying event vector. The 
transverse momentum balance in Z — > ee events is given by 

Ft + Vt + $t = ~Pt c ~ ut- (30) 
Squaring both sides, one finds for the average 

\p e T + Vt + $t I 2 = \VT + ut\ 2 

= n 2 \p T e \ 2 + \u T \ 2 (31) 

assuming again that |pf ec | = /«|p|. e |- The cross term on the right-hand side averages 
to zero since the underlying event vector is randomly distributed with respect to the 
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Z boson recoil system. Figure [12] shows the distribution of + + T | 2 

versus |pf- e | 2 for Z — > ee events. Again, the data demonstrate that there is a linear 
relation between the electromagnetic and hadronic energy scales. The straight line is 
a fit to the data and yields k = 0.83 ± 0.03. This result is consistent with the value 
determined using the other two methods. The intercept of the straight line fit yields 
the magnitude of the underlying event vector \ut\ is 4.3 ± 0.3 GeV. 

Because there was no indication of a non-linear response of the hadronic calorimeter with 
respect to the electromagnetic calorimeter, nor a sign of a measurable offset, the energy scale 
for the hadronic recoil was taken to be strictly proportional to the electromagnetic scale with 
a scale of k = 0.83±0.04, the uncertainty of which was derived from the spread in the results 
among the three different methods. No offset of the hadronic energy scale was included in 
the Monte Carlo model. The effect of a possible non-linearity of the hadronic response was 
included when evaluating the systematic uncertainty on the W boson mass. The only use 
of the EC calorimeter in this analysis was in the determination of the missing transverse 
energy. The hadronic energy scale was determined from the CC-CC (both electrons in the 
CC calorimeter) and CC-EC (one electron in each calorimeter) Z — > ee events. The hadronic 
response was the same as using the CC-derived scale alone within errors which are negligible 
for this measurement. 



C. Resolutions 



1. Electron Energy Resolution 



The electron energy resolution was parameterized according to the relation 



C 2 



s 



+ 



N 



(32) 



'ErJ ' \E. 

where C, S, and N are the coefficients of the constant, sampling, and noise terms, respec- 
tively. The values of the sampling and noise terms were those as derived from test beam 
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data. Smearing in Et rather than in E is used in the sampling term because the resolu- 
tion should become poorer with increasing thickness of the absorber plates at large angles. 
Replacing the usual E with Et compensates for this and allows the coefficient S to remain 
constant over all of the central calorimeter. This observation was confirmed by test beam 
data pB|J31| . The central values utilized in this analysis were obtained from the test beam 



for the central calorimeter and are C = 0.015, S = 0.13 V GeV, and iV = 0.4 GeV. For the 
EC, S = O.lGv/GeV. 

The value for the constant term was determined in situ by fitting the electron energy 
resolution to the observed width of the dielectron invariant mass distribution of the Z — > ee 



events, fixing the width of the Z boson to its measured value of 2.490±0.007 GeV ||19|| . There 
was little sensitivity for small values of the constant term, since for relevant values of Et the 
energy resolution is dominated by the sampling term. A constant term of C = O.OlSitaoii was 
obtained, where the error is statistical only. The uncertainty on the shape of the background 
(discussed below) increased the upper limit on the error to +0.6%. 



2. Electron Angular Resolution 

The polar angle of the electron was determined using the z position of the electromagnetic 
energy cluster in the calorimeter and the z position of the center of gravity of the CDC track. 
The angular resolution used in the Monte Carlo simulation was therefore determined by the 
resolutions on these two quantities. 

The resolution of the calorimeter hit position was determined using electrons from 
W — > ev decays processed through a detailed GEANT Monte Carlo. Because of the read-out 
geometry of the detector, it depended both on the angle of incidence of the electron and its 
cluster z position, z c \ us . It was parametrized as a Gaussian distribution having a width 

Obelus) = (P1+P2 X + (p 3 +P4 X \w\) \z dus \ (33) 

where p x = 0.33183 cm, p 2 = 0.52281 • 10~ 2 cm/degree, p 3 = 0.41968 • KT 3 , and p 4 = 
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0.75496 • 10~ 4 cm/degree. The angle w is the polar angle of the electron (in degrees) as 
measured with respect to the r/=0 axis of the detector. 

The resolution of the z position of the center of gravity of the track was measured from 
Z —>■ ee events using the intersections of the two electron tracks with the beamline. The 
distribution of the difference in z position of the two intersections shows non-Gaussian tails 
which were represented in the Monte Carlo. The simulation generates a resolution on z^ c 
as shown in Fig. 0(a). The model was verified by comparing a Monte Carlo generated 
distribution of the difference in the intersections of the two electrons from Z decays with 
that obtained from the data, and is shown in Fig. |i~3"l(b). 

In the data analysis, the azimuthal angle of the electron was given by the ip angle 
as measured by the CDC. The resolution was taken to be the CDC ip resolution and is 
modeled as a Gaussian distribution with width a ((f) = 0.005 radians. For some Z — > ee 
studies electrons in the end calorimeters were also used. The angular resolutions of these 
electrons were modeled in the Monte Carlo as Gaussian distributions with resolution cr(<p) = 
0.015 radians and a (9) = 0.015 radians. 



The recoil against the vector boson was modeled by assuming it to be a single jet. The 
transverse momentum of the vector boson was smeared using the hadronic energy resolution 
determined both in the test beam and from analysis of jets in situ. It was parameterized as 



3. Hadron Energy Resolution 




(34) 



with resolution parameters C = 0.04, S = 0.80 VGeV and N = 1.5 GeV |2§ . 



D. Efficiencies 



There were two main inefficiencies which affected this measurement: those related to 
the hardware trigger, and those related to electron identification criteria. Both effects can 
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potentially bias the measurement as these particular inefficiencies depend on the kinematics. 
These efficiencies are determined from data as discussed below. 



1. Trigger Efficiencies 

The main data sample was recorded with an on-line filter, which required an electromag- 
netic cluster with > 20 GeV and Ift T > 20 GeV. The trigger efficiency as function of 
the offline electron and missing transverse energies was determined using a single electron 
trigger as well as triggers with lower requirements. After 27% of the running was completed, 
the missing transverse energy calculation in the L2 trigger was changed to use the event 
vertex as measured by the L0 system, rather than the nominal z = value. Therefore, 
two different threshold curves have been used in this data analysis. Both the and $ T 
requirements in the trigger were more than 99% efficient for transverse energies greater than 
30 GeV. 



2. Electron Identification Efficiency 

The recoil of the W boson may affect the electron identification, especially if the recoil 
system is close to the electron. A measure of the event selection biases can be obtained by 
studying identification efficiencies as a function of the quantity u\\, which is the projection 



of the momentum recoiling against the W boson along the electron WA 



u\\ = vT ■ e, (35) 

where e is a unit vector in the electron direction. A bias in the electron identification as 
function of un would distort the lepton spectra. For example, an inefficiency of the 
electron identification at high positive values of u\u when the recoil is close to the electron, 
would result in a softer spectrum. 

The event selection efficiency as a function of u\\ was determined by studying the behavior 



of the energy isolation fraction, /; so , of the electrons in the signal sample. Figure [LJ shows 
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the average isolation versus u\\ for the electrons in the W boson data sample. For negative 
values of un, when the recoil jet is opposite the electron, the isolation is constant. This 
indicates that for these event topologies the recoil system did not affect the electron, as 
expected. For positive values of u\\ the isolation increases linearly with u\\, indicating that 
there was a "halo" of constant energy flow surrounding the direction of the recoil jet. The 
electron identification efficiency was determined by modeling the distribution of the isolation 
variable for different ranges of un as shown in Fig. |15[ The curves are the result of a fit to 
the data using a five parameter functional form. To determine the electron identification 
efficiency as a function of u\\, fits to the isolation distribution were integrated over /; so . The 
fraction of events with / iso > 0.15 constituted a determination of the inefficiency due to the 



recoil jet spoiling the electron signature. The efficiency as function of u\\ is shown in Fig. [16 
where the curve is a parameterized fit. 

The dominant systematic uncertainty stems from the shape of the isolation distribution 
for values of / iso > 0.15, above the trigger restriction. This was addressed by studying W 
boson events in which the electron cluster was rotated in azimuth, re-analyzed, and the 
isolation re-evaluated. The tail of the isolation distribution obtained in this way was well 
described by the fitting function. In addition, when fitting for the isolation distribution 
of the rotated sample, a maximum variation in the efficiency of 1.5% was noted. To be 
conservative, the efficiencies were shifted coherently by two standard deviations of their 
total uncertainties and refit. The band in Fig. [16] shows the resulting uncertainty on the 
efficiency. 



E. Backgrounds 

There was a dual approach to the treatment of backgrounds in this analysis. The process 
W — ► tv — > evvv is indistinguishable from W — > ev and was therefore explicitly included 
in the Monte Carlo event generation. These decays were generated with a 17.9% branching 
fraction for the decay r — > euu, where the r polarization was taken into account. Other 
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backgrounds are characterized by data and were added to the final distributions of the fitted 
variables. The determination of these background contributions is discussed in this section. 

1. Backgrounds toW^ev Events 

The dominant source of background to W — > ev production was standard QCD multi-jet 
production, where one of the jets was misidentified as an electron and there was substantial 
$t from jet energy fluctuations or non-uniform energy response. This background has been 
estimated using data from an inclusive electron trigger that did not impose an isolation 
requirement at the trigger level. The background sample is selected by requiring the same 
kinematic and fiducial cuts as in the W boson event sample but imposing anti-electron iden- 
tification cuts on the EM energy cluster. These anti-electron selections are the combination 
of: 

• /iso > 0.20 

• x 2 > 250 

• (Jtrk > 10. 

For multi-jet background events, it was assumed that the shape of the $t spectrum at low 
$ T was the same independent of the electron quality cuts. The distribution in $ T of the 
background sample was then normalized to the signal sample in the region < $ T < 15 GeV. 
The signal sample was selected from the same trigger by imposing the standard W boson 
selection criteria. The ratio of the number of events with $t > 25 GeV for the signal and 
normalized background distributions was then taken as the amount of background in the 
sample. There was a 0.3% variation in the amount of background due to how the sample is 
normalized and how the background sample was selected. 

As a consistency check, the above procedure was repeated with data taken with an 
inclusive electron trigger that required the EM cluster to be isolated at the trigger level. 
The signal sample was again taken as the events that pass the W boson event selection cuts. 
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The background sample consisted of those events which pass the x 2 > 250 and <Tt r k > 10 cuts. 
Since there was an isolation requirement in the trigger, the background sample does not have 
the anti-isolation cut applied as before. The two methods yielded consistent results. The 
overall background fraction was taken to be (1.6 ± 0.8)%, the average of the two analyses. 
The uncertainty is the total statistical and systematic uncertainty and encompasses the error 
on the two separate measurements. 

Since very few background events survived the kinematic cuts, this method yielded only 
the overall background contribution leaving the shape of the background as a function of the 
transverse mass largely undetermined. Employing the capability of the TRD to distinguish 
electrons, converted photons, and pions a likelihood function was constructed employing the 
energy deposition in the TRD, the track dE/dx in the CDC, and the track cluster match. 
Using an anti-electron criterion based on this likelihood, slightly more background events 
survived the kinematic and acceptance cuts, allowing a determination of the dependence of 
the background as function of the relevant quantities. The data points in Fig. [T7| show the 
calculated transverse mass distribution of the background obtained this way. The line is a 
fourth order polynomial fit. The shape of the background in lepton transverse momentum 
can be described by an exponentially falling spectrum with slope —0.086 ±0.059 GeV -1 and 
—0.129 ± 0.055 GeV -1 for the and Et spectra, respectively. 

The other background that has been considered is the process Z —>■ ee, where one elec- 
tron escapes detection and is not measured (denoted by Z — > e$ ) giving rise to a transverse 



momentum imbalance. This background has been estimated using ISAJET [[37]]. To appro- 
priately model the underlying event in the ISAJET simulation, one minimum bias event was 
vectorially added to the Et for the Monte Carlo data. The overall background contribution 
has been estimated to be (0.43 ± 0.05)%. The M T , p T and E T spectra for this source are 
shown in Fig. [L8|. The E T spectrum does not show a Jacobian edge because the detector 
is hermetic and the energy of the unidentified electron is typically well measured. The 
solid lines for the and p T spectra are from a polynomial fit. The E? spectrum was 
parametrized using an exponentially falling spectrum with slope —0.20 ± 0.03 GeV -1 . The 
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average m for this background is —12.5 ± 0.6 GeV. 



Figure [19] shows the distribution in transverse mass of the dominant background sources 
to the W boson event sample. The background has been normalized to the expected number 
of background events in the data sample. 

2. Backgrounds to Z — > ee Events 

The primary background to Z —>■ ee events came from multi-jet production, with the 
jets fragmenting into a leading tt°. Since the mass is determined from the resonant cross 
section only, a correction also must be made for Drell-Yan and Z7* interference processes. 
These backgrounds were determined as a function of invariant mass and were included at 
the fitting stage. 

The total background contribution was evaluated by fitting the m ee spectrum to a rel- 
ativists Breit-Wigner convoluted with a Gaussian resolution function plus a background 
falling exponentially in m ee . For the mass range of interest, there is no distinction between 
a linear or exponential model of the background. This method yielded a total QCD and 
Drell-Yan background under the Z° peak of 7.4%, with a slope of -0.0447±0.018 (GeV/c 2 ) -1 
for an invariant mass window of 70 to 110 GeV/c 2 . 

The Drell-Yan and Z7* contribution to the total Z boson production cross section was 
determined using an ISAJET simulation. In the mass range 70 < m ee < 110 GeV/c 2 the Drell- 
Yan and Z7* interference terms contributed 3% to the total cross section. The background 
has an exponentially falling spectrum with slope —0.03 (GeV/c 2 ) -1 . The contribution to the 
background from multi-jet sources is thus 4.4%. Both the overall background contribution 
and its shape are in good agreement with the background determination for the cross section 
analysis |38| . 



V. MONTE CARLO SIMULATION 
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A. Introduction 



The W boson and Z boson masses were extracted by comparing measured distributions 
with those generated by a Monte Carlo simulation. To determine the W boson mass the 
relevant distributions are those in transverse mass plus the electron and neutrino transverse 
momenta. For determining the mass of Z bosons the relevent distribution is in the dielectron 
invariant mass. The simulation was accomplished with a generator which produced all of 
the basic processes, incorporated the main features of the D0 detector, and was capable of 
generating tens of millions of simulated events in a few hours. This section describes the 
physics and detector simulation. A comparison between the Monte Carlo simulation and 



the data is presented at the end. Table IV lists the parameters used in the Monte Carlo. 



B. W and Z Boson Production and Decay 

The simulation of W and Z bosons relied on the choice of a model for the physics 
processes involved. This physics was divided into three parts: i) the production model for 
W and Z bosons; ii) the decay of the vector bosons and Hi) backgrounds. For the W 
boson the basic processes generated were W — > ev, W — > tv — > evvv and W — *> 'yeu; for 
the Z boson they were Z — > ee and Z — > ■yee. As discussed in the previous section, all 
backgrounds, except for the W — > tv decay, were not a part of the W or Z simulation and 
were dealt with separately. 



as 



1. Production ofW and Z Bosons 
The triple differential cross section for vector boson production was assumed to factorize 



d 3 a „ da d 2 a 

C— x- — — . (36) 



dprdydm dm dprdy 

Here, C denotes the appropriate constants, y is the rapidity of the vector boson, and the 
products on the right hand side refer to shapes rather than absolutely normalized quantities. 
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The double differential cross section was generated on a grid of pr and y points over the 
region —3.2 < y < 3.2 and < px < 50 GeV/c, in steps of 0.2 in y and 0.5 GeV/c in p T . 
Two choices of the production model, both based on the fully resummed theory of Collins 
and Soper [|3"5[| , were considered. The double differential spectrum as given by Arnold and 
Kauffman ||0| (AK) uses a next-to-leading order calculation for the high p? region |^T| with 



a prescription to match the low and high pt regions. They used fits to Drell-Yan data f42 
which have since been updated. The double differential cross section by Ladinsky and Yuan 
f43fl (LY) employs a different parametrization for the non-perturbative functions describing 
the pt spectrum based on a fit to more recent data. The differential spectra were generated 
for both models using various parton distribution functions as input. Alternative grids 
within the LY model were used, distinguished by a different choice of the non-perturbative 
parameters, (see Appendix A for more details). In order to properly keep track of the 
helicity states for the weak decay, annihilations involving different combinations of valence 
quarks and sea quarks were dealt with separately. The default double differential cross 
section used the LY production model with the MRSA JPJ parton distribution functions. 

After generation of the kinematics of the W boson, the mass dependence of the produc- 
tion cross section was folded in. A relativistic Breit-Wigner line shape was used to model 
the W boson resonance 

da 2 , m 2 

-l—\ m ) TFT- (37) 

dm (m*-M^ 2 + ^f 

where M\y and Tw are the mass and width of the W boson. In pp production, however, the 
mass spectrum differs from the strict Breit-Wigner resonant line shape of the partonic cross 
section due to the variation of parton flux with parton momentum. This mass dependence 
has been calculated by that the differential cross section is given by 



dM 2m f 1 dxi , . / m 2 \ da 



2\ 



= ^(- 2 ) (38) 

m dm 
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with 

Here fq/ p (p)(x) is the probability that a quark or antiquark q in the (anti)proton carries 
a fraction x of the (anti)proton's momentum. In this equation a sum over all qq' pairs 
that lead to W boson production is implicit. The factor JF is a dimensionless quantity 
which will be referred to as the parton luminosity |^5| . It has been parameterized as having 
an exponential mass dependence, e _/3m . The slope parameter (3 has been treated as a 
single number, calculated by evaluating the integral using the available parametrizations 
of the parton distribution functions |46j at a mass of 80 and 91 GeV/c 2 for W boson and 



Z boson production, respectively. The small mass dependence of (3 was included in the 
systematic uncertainty. Table |V| lists the values of (3 which are used for W boson and Z 
boson production for different sets of parton distribution functions. The most recent sets 
which use nearly identical modern input data are the MRSA and CTEQ3M JIT]] sets. The 
relative contributions for vector boson production are listed separately for valence quarks 
and sea-sea quarks. In the event generation the widths of the intermediate vector bosons 



were fixed to their measured values, T(W) = 2.12 GeV || and V(Z) = 2.487 GeV p 



2. Decay of W and Z Bosons 

The W boson decay products were generated in the W boson rest frame with an angular 
distribution depending on which process, valence- valence/sea or sea-sea, is involved. W + 
bosons follow the angular distribution 

d 2 a _ 2 (\doT doV 



~ l-cos0T ■ t:^— + 



dyd cos 9* \2 dy dy J 

(l + cos^.-^- (40) 

where the subscripts v and s refer to valence and sea contributions, respectively, and the +z 
direction is chosen along the proton direction. Here 9* is the center of mass angle between 
the electron direction and the qq axis. 
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The qq — > production cross section at the Z boson resonance is proportional to 

(9 q v+9 q A 2 ) (9v 2 + ^ 2 )(l + cos 2 r) + 

4 cos 0*. (41) 

Because the lepton charge is unmeasured, the cos 0* term averages to zero. The leptons were 
therefore generated with a (1 + cos 2 6*) angular distribution and the uu and dd contributions 
to the production weighted with their respective coupling strengths, g v 2 + g\ 2 . Here g v 
and g\ are the vector and axial-vector coupling strengths of quark q to the Z boson 

g v = Il - 2Q q sin 2 d w (42) 
9 q A = n (43) 

with J3 the third component of the weak isospin and Q q the charge of the quark, if is +2 
for the charge | quarks and -| for the charge -~ quarks. The value sin 2 $ w = 0.2317 |18| 
was used. 

3. Radiative Processes 

Radiative W boson and Z boson decays, q'q — > W — > ez/(7)and qq — > Z — ► ee(7), must 
be properly simulated in the Monte Carlo program to extract correct values for M\y and Mz 
(see Appendix B). Because the ee and eu invariant masses are smaller in these decays than 
the corresponding vector boson masses, the experimentally measured mass distributions 
were shifted toward lower values. 

The rates and distributions in lepton and photon momenta were generated to 0(a) 
following reference ||48|| . Using this calculation, in the decay of the Z boson either of the 
electrons (but not both) may radiate. In W boson decays, the electron or W boson may 
radiate the photon. Approximately 30% ofW^ev decays and 60% of Z — > ee decays had a 
photon of 50 MeV or more in the final state. The calculation does not include processes that 
in the limit of a zero width boson would be considered Wj or Z7 production. Therefore, 
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initial state radiation was not included in the calculation, nor was the production of a virtual 
high mass W boson which decays to an on-shell W boson and a photon. In and Z7 
production, Mw and Mz were correctly obtained from the dilepton invariant masses {ev or 
ee) and the 7 direction was not strongly correlated with that of either lepton. Its presence 
produced a background not fundamentally different from that of other processes. 

In implementing radiative decays in the Monte Carlo simulation, three experimental 
scenarios were considered: i) When the 7 was produced inside the electron cone, taken 
to be a radius of R = 0.2 in r\ — space, the 7 was measured as part of the electron. 
The neutrino momentum, obtained from the missing transverse energy in the event, was 
calculated correctly. Therefore, the invariant mass of the ev system is the W boson mass, 
and the transverse mass and transverse momentum of the W boson was properly calculated. 
ii) If the 7 was far from the electron, that is outside a cone of radius R = 0.4, the photon 
retains its identity. The electron energy was measured correctly, and p\ becomes part of 
the recoil against the W boson, p™ (meas) = pjf (generated) — p}. The transverse mass of 
the ev system was calculated correctly, but was shifted downward because the ev invariant 
mass is smaller than the W boson mass. Therefore, M w was mismeasured. in) If the 7 was 
produced in the region between R = 0.2 and R = 0.4, it alters the shape of the electron 
shower. Isolation and electron identification cuts then resulted in inefficiencies that can 
affect the W boson mass if not properly simulated in the Monte Carlo. 

In the Monte Carlo simulation, the fraction of the electron's energy in the region between 
R = 0.2 and R = 0.4 was generated according to the experimental distribution measured in 
W boson events. The photon energy was added to the electron energy and the event was 
discarded if it failed the isolation cut. If the event survived the isolation cut and the radial 
distance R n between the 7 and e, was less than 0.3, the 7 momentum was added to the 
electron's and the W boson mass was correctly calculated, as in the first case above. If the 
radial distance was greater than R n = 0.3, the 7 momentum was not added to the electron's 
and the reconstructed W boson mass and transverse momentum were shifted downward. 
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C. Detector Simulation 



The production of likelihood templates in M T required large Monte Carlo samples. 
Twenty million generated events were required to sufficiently eliminate effects of statistics in 
the likelihood function. To study the effects of systematic uncertainties many complete anal- 
yses were needed. The combination of these requirements made a fast detector simulation 
essential. 

After production and decay products were boosted into the laboratory frame, the pa- 
rameters whose measurements were described in the previous section were utilized in this 
simulation as follows. 

• The energies of the generated electrons and radiative photons, if they were present and 
retained their identity, were scaled by the measured EM energy scale. The generated 
transverse momenta were then smeared according to the measured resolution, as was 
the generated electron angle. 

• The transverse momentum of the recoil system was taken to be the negative of the 
generated transverse momentum of the W boson, p^ c = —p^. Its magnitude was 
scaled by the product of the measured EM energy scale and the relative response of 
the hadronic and EM calorimeters. Smearing was added according to the jet energy 
resolution. The hadronic content of Z — > ee events and the electrons from the Z boson 
decay were modeled in the same fashion as W — > ev events. 

• The underlying event, denoted by ut{£), was modeled using collider minimum bias 
events, which mimic the debris in the event due to spectator parton interactions and the 
pile-up associated with multiple interactions. The use of minimum bias events properly 
includes any residual energy which might be present from previous crossings as well as 
detector effects. A library of minimum bias triggers was created in bins of luminosity 
in order to correctly simulate overlapping event and noise characteristics of the data. 
Events were chosen according to the distribution of instantaneous luminosities observed 
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during the run as shown in Fig. |2(J (See the discussion in Appendix C.) Figure ^1 
shows the Et and total scalar distributions of the minimum bias events used. The 
average Et is 3.93 GeV with an rms of 2.69 GeV. The mean total scalar Et is 67.1 GeV 
with an rms of 39.8 GeV. (The total scalar Et distribution is shown for completeness 
only, as this quantity is not used in the event modeling.) 

• The generated and smeared recoil hadronic energy vector and the underlying event 
hadronic energy vector were superimposed on one another to form a simulation of the 
total hadronic deposition. 

• The vertex for each generated event was taken to be that of the minimum bias event. 

• The efficiencies and cuts were applied to the smeared quantities. 

1. Underlying Event Discussion 

In the data, the contribution from the underlying event cannot be separated from the 
measured recoil energy. In the simulation of the W events the recoil and the underlying 
event were treated separately. The superposition in the Monte Carlo of the underlying 
event and the production of the W boson and its decay is laced with intricate details. 
Although the average energy deposition per read-out tower in minimum bias events was 
very small, its effect on the W boson mass measurement is of crucial importance mainly 
because the corrections were correlated with the electron direction. Its presence affects 
not only the measurement of the electron energy but also the measurement of the missing 



transverse energy. Equation [H] shows that the neutrino transverse momentum differs from 
the measured missing transverse energy because of the presence of Ut- 

In addition to incorporating the effects of the energy flow of the underlying event in 
the event model, detector effects needed to be taken into account, in particular the effect 
of the zero- suppression. Calorimeter depositions were only read out if the absolute value 
of the magnitude of the energy fell outside the zero-suppression limits. Low energy tails of 
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the electron shower were thus suppressed. A convolution of the shaping electronics and the 
natural radioactivity of uranium caused the pedestal distribution to be asymmetric with a 
long tail towards positive energies. Therefore, even when no particle struck a read-out tower, 
the energy registered for that tower, when read out zero-suppressed, would on average not 
be zero. 

In the following subsection, first the effect of the energy flow of the underlying event 
on the measured energy of the electron will be discussed. The corrections to the measured 
electron energy introduced by the zero-suppression will then be detailed. 

2. Underlying Event Energy and Electron Simulation 

Because the clustering algorithm for electron identification used in this analysis was dy- 
namic, the cluster size can vary from event to event and a description of the underlying event 
contribution to the electron would get rather involved. However, the clustering approach 
was found to be numerically equivalent to a window in (77, <p) space having a constant size 
of 0.5 x 0.5. As this will facilitate the discussion, this window analog consisting of a fixed 
set of 25 towers with (A77 x Aip) = 0.1 x 0.1 will be used to illustrate the size of the effects 
of the underlying event on the electron energy measurement. 

Within the 0.5 x 0.5 window in (77, ip) centered around the electron was contained not 
only the energy of the electron but also the energy from the underlying event. The measured 
electron transverse energy, m p T , is thus given by: 

m F T * Ft + ^ (44) 

where p T refers to the true electron transverse energy, folded with the appropriate resolution, 
and u T 5 is the underlying event contribution inside the 25 towers defining the electron cluster. 
The latter term has been estimated from W events by rotating the electron cluster in azimuth 
and measuring the average energy flow per tower. Care was taken to ensure that the rotated 
cluster was isolated and was not in proximity to any jet activity. The energy flow per 
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tower was found to be 16.8 MeV. The average energy flow under the electron is therefore 
u^? = 25 x 16.8 e = 420 e MeV, with e a unit vector in the electron direction. 

This contribution has also been determined from minimum bias events, spanning an 
appropriate range in luminosity. An average energy flow of 15.3 MeV per tower was found. 
The difference of 1.5 MeV between the two methods is attributed to the presence of the W 
recoil. A value of 16.8 MeV per tower has been used in the simulation. The uncertainty on 
the average energy flow is reflected in the systematic uncertainty due to this source. 

To each of the two terms in Eq. fl44]) , a correction needed to be applied due to zero- 
suppression. Under normal running conditions, calorimeter cells were not read out if the 
signal was within 2cr of the mean pedestal for that channel; that is, the read-out was zero- 
suppressed. As a consequence, the tails of the electron shower which fell within the zero- 
suppression limits Uj> were lost. The average energy that was lost below the 2a zero- 
suppression threshold was estimated to be u^ zs = —152 e MeV for electrons from W decays, 
using a detailed GEANT simulation. 

Because the absorber medium in the calorimeter is uranium, which is a natural /^-emitter, 
the pedestal distributions were asymmetric. Additionally, some asymmetry in the pedestal 
distributions was introduced due to the shaping electronics P^ ]. Therefore, even when no 
particle strikes a read-out tower, the zero-supressed energy read out for that tower was 
on average not zero. This zero- suppression contribution has been studied by analyzing non 
zero-suppressed minimum bias events. By comparing the energy per 0.1 x 0.1 read-out tower 
measured in these events to the energy that results after applying the zero suppression offline, 
the energy per read-out tower of EM and the first FH layers, was 7.55 MeV higher than 
in non zero-suppressed events. It should thus be realized that the average energy flow of 
16.8 MeV per tower, derived above, has two contributions. The first contribution is from the 
true energy flow in the event, determined to be 9.23 MeV per tower. The second contribution 
is an artifact of the zero-suppression, due to the asymmetric pedestal distributions, which 
adds an energy of 7.55 MeV per tower to the read-out. 

As mentioned before, a minimum bias event was used to model the event underlying 
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the W boson. The presence of the electron from the W decay affected the energy flow in 
the underlying event. Notably, the read-out towers occupied by the electron had a very 
large energy deposition and therefore were not affected by the zero-suppression correction. 
For the W data, the electron occupied on average 8 ± 3 towers. Therefore, applying the 
zero-suppression correction to all 25 read-out towers of the electron cluster, which has been 
assumed above, is incorrect. This was corrected by applying the correction to only the 
17 channels within the cluster that on average were zero-suppressed or, equivalently, by 
subtracting out the zero-suppressed pedestal energy from the 8 cells that on average were 
read-out with the electron. Thus, a correction Uj? s = — 8 x 7.55 e MeV needed to be applied 
to the energy flow under the electron, w| 5 . 

To summarize, the measured electron transverse energy, in MeV, is given by 

= Vt - 152 e + 25 x (9.23 + 7.55) e - 

8x7.55e (45) 

with 

• p|, the true electron transverse energy folded with its resolution; 

• Uj< w the energy of the tails of the electron shower lost due to the zero-suppression, 
determined to be — 152 e MeV; 

• Ut 5 the energy flow from the underlying event within the 0.5 x 0.5 window in (77, <p) 
defining the electron cluster, given by 25 x (9.23 + 7.55) e MeV; and 

• the correction to the energy flow of the underlying event due to the presence of 
the electron which corrects for the zero-suppression effect of the underlying event for 
the towers occupied by the electron, —8 x 7.55 e. 

When all of these effects were taken into account, an addition of an average of 207 MeV 
to the generated electron along the electron direction was required in order to correctly 
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simulate the measured electron px- 



3. Underlying Event Energy and Recoil Energy 

The measured recoil energy in the detector is a combination of the true recoil of the W 
boson and the contribution of the underlying event. In the simulation the true recoil of the 
W boson was taken to be p T ec and the underlying event was simulated using a minimum 
bias event. Therefore the measured recoil was given by 

J™ ~ -pf + u T . (46) 

The underlying event vector ut was taken to be the sum of the Et of all calorimeter cells 
in the minimum bias event. However, a correction needed to be applied to the underlying 
event energy vector due to the presence of the electron in W events. Recall that in the data 
analysis the recoil momentum was determined by subtracting the electron transverse energy 
from the total measured transverse energy in the event. Therefore the energy flow under the 
electron from the underlying event, pointing along the electron, should be subtracted from 
Ut- In the simulation, the recoil was thus calculated as 

->rec _ , 7 725 

mPr — Pt + u t — u t 

= - ffl + u T - 25 x (9.23 + 7.55) e . (47) 

Note the absence of the Ut term, which does not need to be applied here since ut is from 
a minimum bias event in which no high pt electrons are present. 

The underlying event model and the resolution in p T ec has been verified using the 77 
imbalance in Z boson events, defined previously. Since the magnitude of the Et in minimum 
bias events was of the same order as that of the pt of the vector boson, the width of the 
distribution of the 77 imbalance (see Fig. pH|(b)) was very sensitive to the underlying event 
contribution. The rms of the 77 imbalance distribution in Fig. |TU|(b), after the correction 
for the hadronic energy scale has been applied, is a = 4.44 ± 0.18 GeV. This is the band 
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shown in Fig. By varying the number of minimum bias events in the Monte Carlo 
that mimic the underlying event, the width of the rj imbalance determined the number of 



minimum bias events to be added in the simulation. The points in Fig. show the Monte 
Carlo predicted widths as function of the number of minimum bias events. The number of 
minimum bias events preferred by the data was iV min . bias = 0.98 ± 0.06 events. Since this 
number is consistent with 1.0, one minimum bias event was used to model the underlying 
event in W and Z boson production. 



4- Underlying Event and the Neutrino Momentum 

The neutrino momentum is a derived quantity which follows directly from the electron 
and recoil measurements: 

$ T = -JT - Jt 

= -Pt c - Pt -ut- uu s (48) 

where ut zb = + «ji ■ Note that ut — represents the energy vector of the underlying 
event with the region that the electron occupies excised. 

There are two equivalent ways to view the effect of the underlying event. If one uses for 
the neutrino momentum the second line above, then the measured electron energy, including 
the contribution from zero suppression and the energy from the underlying event, appears in 
the neutrino and the electron in W decays and in both electrons in Z decays. This correction 
then cancels in the ratio of the two masses. Then what is important is the amount of the 
underlying event energy which should be excluded from the determination of the W boson 
recoil energy because it is inside the electron cluster. Alternatively, if one examines the 
expression for the neutrino momentum given in the third line above, only the total recoil 
momentum and the total underlying energy enter. The zero suppression correction is still 
irrelevant, appearing in the neutrino, the W electron, and the two electrons from the Z boson 
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decay. Now the correction to the electron energy from the energy flow from the underlying 
event that appears inside the electron cluster does not cancel completely in the M w /M z 
ratio. 

The missing transverse momentum differs from the neutrino momentum because of the 
presence of u T - This effect has no counterpart in Z boson decays and it changes the measured 
transverse mass and must be properly modeled. As described above, in the Monte Carlo 
simulation ut was obtained from minimum bias events. If there were a biased region of 
the calorimeter which made ut directional, this effect would be accounted for in the Monte 
Carlo events. Although the above is dependent on properly extracting small energies in the 
calorimeter, many of the effects cancel in the ratio M w /M z . 

D. Application of efficiencies 

After simulating the vector boson event kinematics, the efficiencies of the trigger as well 
as the electron identification efficiency as a function of u\\ were applied, using the measured 
kinematic quantities. Fiducial cuts in 77 and <fi were made as in the data. Using the measured 
quantities, the transverse mass was calculated and the same selection criteria as in the data 
were applied: m T > 50 GeV/c 2 ; E e T > 25 GeV/c; $ T > 25 GeV; and Vt < 30 GeV/c. 

E. Comparison of Data with Monte Carlo 

Comparisons of various distributions of the simulated quantities with data are shown 
in this section. The distributions comparing the data and the results of the simulation are 
area normalized. The Monte Carlo was generated at the final W boson mass value of this 
analysis obtained from the transverse mass fit. In the comparisons the data are generally 
shown as points with statistical errors; the simulation is shown as the histogram. 
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1. Characterization of the W — > ev Candidates 

The primary measurables in W — > ev events are the energy and direction of the electron 
E and the transverse momentum of the recoil p^ c . In addition, there are a variety of derived 
quantities which are especially sensitive to the presence of inefficiencies or biases which serve 
as important checks. The comparison between the data and the Monte Carlo simulation 
for W —>■ ev events in the electron polar angle cos(# e ) and the transverse momentum of 
the W boson p^ are shown in Figs. and There is reasonable agreement between the 
simulation and the data in both distributions. 

Because of its strong correlation with the lepton transverse momenta, un, defined pre- 
viously, is an important quantity. As was noted in Section IV D2, a bias in u\\ distorts 
the available momentum phase space of the leptons and results in a softer or harder lepton 
Pt spectrum, depending on that bias. Since un involves both the electron identification 
efficiency and the hadronic energy scale, it is advantageous to study the distribution in 
the angle between the recoil system and the electron, as well as a distribution in un itself. 



Figure |25] shows the distribution in ip e i — ip rec . 

Note that for small , assuming perfect electron identification, the W boson recoil would 
be distributed uniformly in (p around the electron direction. However, the distribution in 
fei — free is asymmetric. There are two sources for this asymmetry. The dominant effect is 
simply the kinematics of W — *> ev decays. For transversely boosted W bosons, on average 
the electron carries away p e T m pj! '/2 along the p^f direction, having a magnitude of ~ My/ jl 
for small values of pYp . This implies that (u\\) ~ —{p^ 2 )/Mw- Since the mean value of p? 



is approximately 9 GeV/c (see Fig. |24j) , (wii) is about —1 GeV and the distribution in the 
difference in azimuthal angle of the electron and the recoil tends to favor negative values 
of un. The second effect which enhances the asymmetry is due to a decrease in electron 
identification efficiency as function of u\\. The value of u\\ is an indication of the proximity 
of the recoil jet to the electron. For high positive values of u\\ the recoil jet is close to 
the electron and can spoil its signature. The observed excellent agreement between the 
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simulation and the data indicates that the event kinematics and the electron identification 
efficiency are modeled adequately. 

Figures [26] and show the correlation between (u\\) and p^ and p^. An important 
feature of the transverse mass is that, unlike pf> and pj<, Mt is relatively uncorrelated with 
u\\ as shown in Fig. EM This shows clearly one of the advantages of using the transverse 
mass to obtain the W boson mass. The correlation between u\\ and p^ is shown in Fig. 29. 



Figure |30| shows the distribution in un itself. Note that there has not been a subtraction 
for background. The mean value of u\\ for the data is (u\\) = —1.19 ± 0.08 GeV whereas the 
simulation gives (u\\) = —1.13 ± 0.02 GeV. An average correction for the QCD and Z — > ee 
background has been applied to the value just quoted for (uii) for the Monte Carlo. 



The distribution of u±, Fig. |31], defined as the projection of the recoil jet onto the axis 
perpendicular to the electron direction, is a measure of the resolution of the recoil system. 
Its mean value is close to zero, as expected. For the data (u±) = 0.025 ± 0.087 GeV with 
an rms of 7.4 GeV; the simulation gives (u±) = 0.024 GeV with an rms of 7.5 GeV. 



2. Characterization of the Z — * ee Candidates 

The measured quantities in Z — > ee events are the energy and direction of both electrons 
and the transverse momentum of the recoil system. Equally important are the determination 



of derived quantities of the Z boson kinematics. Figures 32 show the comparison in 



electron energy E e i and the transverse momentum distribution from the recoil system, p r rf c . 



VI. FITTING PROCEDURE 

The Monte Carlo event generation was performed for 21 equidistant mass values binned 
at intervals of 100 MeV/c 2 in the transverse mass for W boson spectra, 200 MeV/c 2 in 
invariant mass for Z boson spectra, and 100 MeV/c for the transverse momentum spectra. 

An unbinned maximum likelihood fit was used to determine the vector boson mass using 
the normalized Monte Carlo spectra as templates. The log-likelihood was calculated for the 
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data for the 21 different generated masses. Since the templates were binned whereas the 
data were unbinned, a quadratic interpolation between adjacent bins in the templates was 
performed. The log-likelihood values for the 21 different vector boson masses were fit to a 
parabola and the minimum was taken to be the fitted mass value. A decrease of half a unit 
in the log-likelihood is the quoted single standard deviation statistical uncertainty. 

The likelihood distribution need not be Gaussian, depending on the range of the param- 
eter fit, the intrinsic shape of the spectrum and the resolution function. This is particularly 
true for spectra with a sharp edge like the Jacobian peak in the distributions considered 
here. Both quadratic and cubic polynomial fits were performed to the log-likelihood. The 
differences were small and for all results presented here, the values from the quadratic fit 
are quoted. 

Any Monte Carlo-based fitting procedure should satisfy the requirements that, if the 
procedure is applied to an ensemble of Monte Carlo generated data samples, it returns the 
input values with which the events were generated and, secondly, that the rms spread of the 
values for the fitted parameter be consistent with the mean statistical uncertainty of the 
fit to each individual data sample. This was done for an ensemble of 125 generated data 
samples of 8000 events each. The average statistical error for each of the three different W 
boson mass fits is: 8(M T ) = 130, 5(p e T ) = 183 and 5(p v T ) = 248 MeV/c 2 , respectively. The 
average fitted mass values are M W (M T ) = 80.410 ± 0.013, M w (p e T ) = 80.398 ± 0.017 and 
M w [f T ) = 80.420±0.021 GeV/c 2 , in good agreement with the input value of 80.400 GeV/c 2 
within the statistical accuracy of the generated templates. They are consistent with the rms 
spread of the distribution of the fitted masses, rms(Mr) = 145 ± 9, rms(p^) = 188 ± 12 and 
rms(p^) = 237 ± 15 MeV/c 2 , respectively. Figure [M] shows the distribution of fitted mass 
values and fit uncertainty for W bosons as obtained from a fit to the transverse mass for 
this ensemble. 

As discussed in the previous section, backgrounds were not included in the event simu- 
lation. Their effect on the mass determination was taken into account through inclusion of 
the shape of the background spectrum in the likelihood distributions. The background was 
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properly normalized to the expected background fraction in the relevant fitting range. All 
results were corrected for backgrounds. 



A. Results of Z Boson Mass Fits 

The dielectron invariant mass spectrum for the central-central (CC-CC) event topology, 
with the corresponding best fit of the templates to the data, is shown in Fig. |35|. The events 
in the mass range 70 < m ee < 110 GeV/c 2 were used to extract the Z boson mass. The final 
measured Z boson mass for events which require both electrons in the central calorimeter 
is: 

M z = 91.070 ± 0.170 GeV/c 2 . (49) 



The error is statistical only. Figure 35 also shows the relative likelihood distribution and 



signed a/x 2 " of the fit for central-central electrons. 

B. Results of W Boson Mass Fits 

The W boson mass was obtained from fits to the transverse mass of the W boson, Mt 
(Fig. |36|), the electron px (Fig. ^) and the neutrino px spectrum (Fig. |38[). The transverse 



mass fit was performed over the range 60 < M T < 90 GeV/c 2 , which contains 5982 events. 
Placing the lower edge at 60 GeV/c 2 removed most of the QCD background. Since the 
probability for finding events in the very high transverse mass tail was small, relatively 
small fluctuations in the number of observed high transverse mass events can significantly 
affect the fitted mass. Given that the high transverse mass tail of the QCD background 
was rather poorly known, a high Mt cut of 90 GeV/c 2 was also imposed. A transverse 
momentum range of 30 to 45 GeV/c was used for fits to the transverse momentum spectra. 
There were 5520 events in the fitting range for the electron transverse momentum spectrum 
and 5457 events for the neutrino transverse momentum spectrum. It should be noted that 
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the fitting windows were placed on "uncorrected" energies, that is, electron energies which 
had not been scaled as described in Section IV. 

The final fitted masses from the three spectra are 

80.350 ± 0.140 GeV/c 2 (50) 
80.300 ± 0.190 GeV/c 2 (51) 
80.045 ± 0.260 GeV/c 2 (52) 

The errors are again statistical only. Note that the W boson mass determination using 
the transverse mass is the most precise. After taking into account the small offset, which 
resulted in a change of the W boson mass of 5 MeV/c 2 as described in Section IV A, the 
measured mass ratio is 

M w /M z = 0.88114 ± 0.00154 (53) 
where the error is statistical only. 

VII. SYSTEMATIC SHIFTS AND UNCERTAINTIES 

In this analysis, the W boson mass was obtained from a fit to the spectrum in transverse 
mass defined in Eq. ([12]). The Z boson mass was obtained from a fit to the spectrum in 
invariant mass of the two electrons, defined in Eq. (|H]). In this section the uncertainties in 
the measured masses that could arise from mismeasurements of the terms in these equations 
are described. Note that the errors quoted will be those for the measured W boson mass 
which is extracted from the ratio of the fitted W boson and Z boson masses; correlations 
between the two masses have been taken into account. 

Unless otherwise noted, the determinations of the shifts in mass due to the various un- 
certainties have been obtained through Monte Carlo studies and are labeled "Monte Carlo" 
in the tables. In these studies, high statistics Monte Carlo event samples were generated 
with the parameter in question varied within its allowed range. These samples were then fit 
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M W (M T ) = 
M w (p e T ) = 
M W (P») = 



to the templates with the nominal settings to determine the systematic error. The errors on 
these shifts reflect the statistical error on the simulation. The sensitivity, 9I ^p , where P is 
the parameter that has been varied, was determined from a linear fit to the shifts in mass 
over a representative range around the nominal value of the parameter. Values in the tables 
labeled "Data" are the shifts in mass when the data are fit to a template in which one of 
the parameters deviates from its preferred value, with the others unchanged. No error is 
quoted for these data shifts, since it would be meaningless. 



A. Electron Energy Scale Uncertainty 



As discussed in Section [IV A| , many systematic effects due to the calorimeter scale which 
are common to the measurement of both the W and Z bosons cancel in the ratio of their 
masses. However, there are small effects that can bias the measured Z boson mass in ways 
which do not cancel in the ratio, Mw/Mz, and they are discussed in the next section. 



1. Uncertainties in Mz 

The first source of a possible bias in the Z boson mass measurement is the background 
under the Z boson resonance. The nominal multi-jet background in the Z — > ee sam- 
ple and the Drell-Yan contribution caused a shift in the reconstructed Z boson mass of 
+39 ± 12 MeV/c 2 . The uncertainty on this correction has been estimated by varying the 
slope of the background which resulted in a change in the overall background level from 3.2% 
to 8.2%. Such a variation in the background results in a variation of 20 MeV/c 2 in M z , 
which was taken to be the systematic uncertainty on the Z boson mass from the background 
contribution. Other uncertainties arose due to parton distribution functions, radiative cor- 
rections, and a small fitting error. Among these, the change in parton luminosity for the 
different parton distribution functions was most significant. Varying the parton luminosity 
slope (3 within the range given by the various parton distribution functions considered in 
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this analysis, 1.030 x 10 2 < (3 < 1.113 x 10 2 , along with the other effects results in an 
overall 35 MeV/c 2 uncertainty in the Z boson mass. 



2. Total Mw Uncertainty Due to Electron Scale 
As was noted in Section [IV A|, the largest contribution to the overall scale uncertainty 



was due to the number of Z boson events. This statistical component was 150 MeV/c 2 . 
In addition, the uncertainty due to the possible nonlinearity in the calorimeter response as 
determined by the combined rriM, m T o, and Mz analysis (related to the uncertainty in 5) 
was assigned as 25 MeV/ c 2 . Combining these in quadrature with the systematic uncertainties 
just discussed resulted in the overall scale uncertainty assignment which is rounded up to 
160 MeV/c 2 . 



B. Uniformity of Electron Energy Response Uncertainty 

The data were corrected for the observed azimuthal variations in energy response of the 
different calorimeter modules, reducing the error from this source to a negligible level. Any 
residual non-uniformity in response was taken into account through the constant term in 
the energy resolution. 

A non-uniform response in rj, however, can introduce a bias in the mass determination, 
arising from the fact that the kinematic distribution of electrons from Z boson decays differs 
from that in W boson decays. The electrons from Z boson decays have a different average rj 
than the electrons from W boson decays, even when event samples are very large. Moreover, 
a non-uniformity can distort the differential distributions. A possible r\ dependence of the 
calorimeter response will thus not cancel in the ratio of the two masses. 

To address this, the response of the different rj regions of the detector were scaled in the 
Monte Carlo with respect to the nominal uniform response. Two sets of scale factors were 
used, corresponding to the response of two EM modules measured in the 1991 test beam. 
These scale factors were applied in discrete steps in 77, following the read-out geometry of 
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the calorimeter, and varied from 0.985 to 1.013 over the central pseudorapidity range. The 
observed shifts in fitted mass are listed in Table [VT]. Assuming the rj response of the test 
beam modules typified the variation in uniformity, a systematic uncertainty on the W boson 
mass from the transverse mass fit of 10 MeV/c 2 was assigned due to this uncertainty. 



C. Electron Energy Resolution Uncertainty 

The electron energy resolution in the central calorimeter was parameterized as discussed 
in Section IV C3. Most effects which degrade the resolution affected the resolution function 
constant term. For example, spatial non-uniformities in the detector response and electronics 
gain variations contributed to the constant. The sampling term varies very little, from 1.9- 
2.4%, as the electron p? is varied over the range 30-45 GeV/c. Therefore changing only the 
constant term and noting the change in the W boson mass was sufficient to accommodate 
most sources of uncertainty in the energy resolution. 

To study the dependence of the W boson mass on the resolution, the constant term was 
varied in the Monte Carlo simulation. The W boson mass increases if a resolution smaller 
than actually exists in the data is used in the Monte Carlo. Better resolution in the Monte 
Carlo results in a sharper Jacobian edge and the fitted mass shifts upward to accommodate 
the larger resolution tail in the data. The transverse mass distribution was most sensitive, 
since the Jacobian edge was best preserved. For the px spectra the edge is smeared, due in 
part to the transverse boost of the W boson. Table [VI I| lists the changes in W boson mass 



for all three fits when varying the constant term by 0.5% from its nominal value of 1.5%. 
An uncertainty in the measured W boson mass for the transverse mass fit of 70 MeV/c 2 was 
assigned according to this variation. 



D. Electron Angle Uncertainty 



The electron polar angle is defined by the position of the electromagnetic cluster in the 
calorimeter and the position of the cog of the CDC track. Recall from Eq. (|18|) in Section 
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Ill that a scale factor, ocdc, was applied during the data analysis to correct the bias in the 
z position of the cog of the CDC track. The uncertainty in the W boson mass due to the 
uncertainty in q;cdc has been determined by applying varying scale factors to the z position 
of the CDC cog in the W boson and Z boson data and fitting to the standard templates. By 
varying the CDC scale factor around the nominal value within its tolerance of 0.002 for the 
W boson and Z boson data sample simultaneously, the uncertainty on the W boson mass 
was determined to be 50 MeV/c 2 . 



E. Hadronic Energy Scale Uncertainty 

The energy scale of the vectors «y and p? , which both include hadronic energy, was not 
the same as the scale of p^, which contains only electromagnetic energy and was calibrated 
by the Z boson mass. The relative hadronic to electromagnetic energy scale is set using Z 
boson events and the scale obtained is k = 0.83 ± 0.04 . The sensitivity of the measured 
W boson mass was obtained by varying the value of « within its uncertainty in the Monte 
Carlo generation of the templates. The 0.04 variation in hadronic energy scale produced a 
50 MeV/c 2 uncertainty on the W boson mass from the transverse mass fit, where an increase 



in the scale factor resulted in an increase of the measured W boson mass. Table |VIII| lists 
the change in W boson mass when varying the hadronic energy scale factor by 0.04 from its 
nominal value for all three fits. The mass obtained from the p^ fit was affected by a change 
in the hadron energy scale through the electron identification efficiency as function of un. 



F. Hadron Energy Resolution Uncertainty 

The resolution in p^ had two components: the energy resolution of the recoil jet which 



is aligned with the recoil direction P5 |, and the underlying event vector Ur which was 
randomly oriented with respect to the recoil. In the Monte Carlo the recoil momentum 
p^ c was simulated by assuming it is a jet with resolution ahad/E = 80% /VE as discussed 
above. All of the uncertainty due to this quantity was presumed to be accounted for through 
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variations in the sampling term alone. The second contribution, that from ur, dominated 
the overall resolution in p^f . It was obtained directly from the experiment using minimum 
bias events chosen at the proper luminosity to simulate the underlying event. 

The data constrained the number of minimum bias events to A^ min .bi as = 0.98 ± 0.06. 
The nominal value used in the simulation was 1.0. The change in W boson mass for various 



values of the number of underlying events is listed in Table [LX|. This includes the effect 
of resolution broadening and the neutrino scale shift which results from changing ut- The 
application of the randomly oriented underlying event has the effect of adding an azimuthally 
symmetric component to the overall resolution for the total hadron energy vector. The 
systematic uncertainty on the measured W boson mass due to the uncertainty on the number 
of minimum bias events is 60 MeV/c 2 for the transverse mass fit. 

The mass determined from the p^ spectrum was, within errors, not affected by the hadron 
energy resolution. The W boson mass determined from the other two spectra would increase 
if a smaller average number of minimum bias events underlying the W boson were used in 
the Monte Carlo since the resolution improves. 

The jet energy resolution also contributed to the uncertainty attributed to the overall 
hadronic energy resolution. Varying the sampling term in the jet energy resolution from 0.6 
to 1.0 changes the W boson mass by 65 MeV/c 2 , which was taken to be the systematic error 
due to this source. Table |X| lists the change in the mass from the different fits when varying 
the sampling term of the hadronic energy resolution. 

G. Energy Under the Electron Uncertainty 

The measured electron energy not only consisted of the electron energy itself, smeared 
by the detector resolution, but also included a contribution from the underlying event. In 
addition, there was a bias in the electron energy due to zero-suppression in the readout 
electronics. Following the discussion in Section |V C| the measured electron px was modeled 



as a combination of four terms, p%, Uj, , uj ; , and . The additional contributions to 
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the electron energy point, to a good approximation, along the electron direction with the 
magnitude of 207 MeV. The uncertainty on this has been estimated to be approximately 
50 MeV. The measured neutrino momentum can be written in two equivalent ways: 

$T = ~Pt C ~ m PT - [UT - U^} 

= -p T T ec - p° - u T - u Tzs (54) 

Using the second equation, the total recoil momentum and the total underlying energy enter 
in the calculation of the neutrino momentum. Both were well determined by the W boson 
and Z boson data. Using this approach, the overall uncertainty derived from the measured 
electron energy in a manner which did not completely cancel in the ratio Myy/Mz- The zero 
suppression correction here was quite small, since it contributed to the neutrino and the W 
boson electron, as well as the two Z boson electrons. 

Using the first equation the measured electron energy appeared in both the measured 
neutrino momentum and the measured electron momentum for W boson decays and in both 
electrons for Z boson decays. The correction to the electron energy then canceled completely 
in the ratio of the W boson and Z boson masses. What is important is how much of the 
underlying energy and W boson recoil energy should be excluded from the event for the pj? c 
determination, because it was inside the electron cone. The method used to determine the 
uncertainty on the W boson mass from the contribution due to energy under the electron 
followed this approach. 



Three effects were identified that contribute to this uncertainty. Figure |39| shows the av- 
erage transverse energy flow in an EM tower plus the first FH layer versus tower index (i'L). It 
is seen that the energy flow was constant in rj within 0.5 MeV for the central calorimeter. In 
the Monte Carlo a uniform distribution was assumed and the deviation of a flat distribu- 
tion from that shaped like the data contributed an uncertainty of approximately 20 MeV/c 2 
on the W boson mass. 

The second source of uncertainty stems from the fact that the underlying energy in W 
boson events was measured to be 16.8 MeV per tower in the EM plus FH1 layers, whereas 
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minimum bias events yielded 15.3 MeV. In the Monte Carlo an energy flow of 16.8 MeV was 
assumed. This difference of 1.5 MeV, most likely due to the presence of the W boson recoil, 
was treated as an uncertainty on the mass which is equal to (25 x 1.5)/2 ~ 20 MeV/c 2 . 

The third source is due to the uncertainty on the number of towers to be excluded from 
the Et of the underlying event. In the Monte Carlo, a region of 5 x 5 = 25 towers was 
excluded. In the data, the number of towers used by the electron in the clustering algorithm 
varied event by event. This uncertainty on the W boson mass was evaluated by repeating 
the analysis using another electron clustering algorithm that removed this error completely 
(see Section [VIII B ) . The difference in W boson mass between the two electron clustering 



approaches led to a 20 MeV/c 2 uncertainty due to this effect. These three uncertainties were 
summed in quadrature to obtain the total uncertainty on the W boson mass of 35 MeV/c 2 
due to the uncertainty in the energy flow underlying the electron. 

H. Production Model Uncertainty 

In the generation of the W boson and Z boson events a theoretical model for the vector 
boson transverse momentum and rapidity spectrum was used. This production model had 
an uncertainty associated with it which led to an uncertainty in the measured W boson 
mass. Since parton distributions and the spectrum in are correlated, this correlation 
was addressed in the determination of its uncertainty on the W boson mass. To constrain 
the production model, both the measured spectrum as well as the published CDF W 



boson charge asymmetry data |50| were used. 

The parton distribution functions were constrained by the CDF measured W boson 
charge asymmetry data. To accommodate the variation allowed by the asymmetry data while 
at the same time utilizing the available data from all other experiments, new parametriza- 
tions of the CTEQ3M parton distribution function were obtained The fit used to 



obtain these parametrizations included the CDF W boson asymmetry data with all data 
points moved coherently up or down by one standard deviation. These parametrizations 
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will be referred to in the following as "asymmetry high" and "asymmetry low", respec- 
tively. Figure [ID] shows the relative change in the theoretical pf spectrum for these new 
parametrizations of the CTEQ3M parton distribution function with respect to the nominal 
spectrum. 

The pt spectra of the vector bosons were most sensitive to variations in the parameter g 2 , 
which describes the Q 2 dependence of the parametrization of the non-perturbative functions 
(see Appendix |Aj). Figure |4l| a shows the change in the pf spectrum when the parameter g 2 
is varied significantly from its nominal value. Note that for low pt, the cross section varies 
by approximately a factor of two. Figure [4l] b shows the constraint on g 2 by the Z boson 
data as given by a simple x 2 test. For the estimate of the uncertainty on the W boson mass, 
the range for g 2 was limited to —2a < g 2 < 4cr, which are conservative bounds in agreement 
with the Z boson data. 

To assess the uncertainty due to parton distribution functions and pYjt input spectrum, 
the change in W boson mass was noted when varying both the parton distribution function, 
as determined by varying the measured W boson charge asymmetry, and the g 2 parameter 
simultaneously. The results of the change in W boson mass are listed in Table KIl A total 
error on the W boson mass of 65 MeV/c 2 has been assigned due to the uncertainty on the 
parton distribution functions and the input p^ spectrum. 

The change in W boson mass obtained from high statistics Monte Carlo studies for 
different parton distribution functions, compared to the nominal MRSA parton distribution 
function is shown in Table |XII| . An uncertainty of 50 MeV/c 2 in the measured W boson 
mass could be attributed to the choice of parton distribution function. Note that this 
uncertainty is only listed for completeness. The more conservative estimate, varying both 
the parton distribution functions and the p^ spectrum simultaneously, was taken as the 
final uncertainty due to these sources. 

Finally, the finite width of the W boson was taken as 2.1 ±0.1 GeV and the effect on 
the W boson mass due to its uncertainty was found to be 20 MeV/c 2 . 
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I. Background Uncertainty 



The presence of background caused a bias in the determination of the mass. The shift 
in mass has been determined by including the nominal background spectra in the likelihood 
templates. Systematic uncertainties arose due to the uncertainty on the overall background 
contribution and the shape of the background spectrum. 

The QCD multi-jet background contribution to the signal sample is (1.6 ± 0.8)%. The 
contribution of Z — > ee, in which one electron is not identified, is (0.43 ± 0.05)%. The 
presence of these backgrounds introduced a shift in measured mass of +33 MeV/c 2 and 
+4 MeV/c 2 , respectively, for the transverse mass fit. The background levels have been 
varied within the quoted uncertainties. The shape of the QCD multi-jet background for 
the transverse mass distribution was varied as shown by the curves labeled "excursions" in 
Fig. |42|. Similarly, the shape of the Z — > ee background was varied. The total systematic 
uncertainty on My/ due to the variations in the QCD and Z — > ee background is 30 MeV/c 2 
and 15 MeV/c 2 , respectively. An overall uncertainty of 35 MeV/c 2 has been assigned to the 
uncertainty in the background. 

J. Radiative Decay Uncertainty 

The parameters used in the modeling of radiative decays were the minimum separation 
between the electron and photon for the photon to retain its identity R ei and the minimum 
energy of the radiated photon _E™ m . The uncertainty in the value of these parameters 
led to an uncertainty in the measured W boson mass. Uncertainties can also arise from 
inefficiencies caused by the photon affecting the electron shower shape, the effect of upstream 
material on the energy measurement of photons and from theoretical uncertainties. 

The electron photon separation parameter R ei was varied by ±0.1 from its nominal value 
of 0.3 and the effect on the W boson mass was noted. From this, an uncertainty of 10 MeV 
on My/ was determined. In a second independent analysis the correlation between the effect 



66 



of a photon on the isolation as well as the topological requirements was taken into account 
through a full detector simulation. The four-vectors of the decay products from radiative 
decay events were input to the GEANT simulation. The events, processed using the standard 
reconstruction algorithms, were then subjected to the same selection criteria as the data and 
electron identification efficiencies were determined as a function of E% and R n . Modeling 
the resulting variation of the efficiencies determined in this fashion in the Monte Carlo led, 
again, to an uncertainty on the W boson mass of 10 MeV/ c 2 which is the same as that found 
in the other method. 

The dependence of My/ on E™ m was negligible. The choice of E™ m = 50 MeV was suffi- 
ciently low that within the accuracy of the measurement it was insensitive to this parameter. 

In the modeling of radiative decays, only order a EM corrections to the lowest order dia- 
grams have been considered and processes in which two or more photons are radiated have 
been ignored. Also, initial state radiation and finite lepton masses were not included in the 
calculation. This effect has been estimated to be 10 MeV/c 2 and confirmed by a recent 
theoretical calculation IE2] . 



Since the effect of radiative decays was large and changed the W boson and Z boson 
masses in a way that did not cancel in the ratio, it was important to also evaluate the effect 
when the photon is produced by bremsstrahlung in the central detector. For the photon to 
have an effect on the measured W boson mass, it must be separated from the electron in 
(77, <f) space by at least R = 0.2. 

The probability for radiating a photon is very strongly peaked at small angles (see Ap- 
pendix B) , with very little dependence on the fraction of the electron's energy carried by the 
photon ||53|| . The photon never separates from the electron beyond a cone of 0.2 by radiation 
alone and therefore external bremsstrahlung has no effect on the W boson mass. 

As shown in Appendix B, the electron and photon can also separate if the electron 
undergoes multiple scattering through a large angle. The angles resulting from multiple 
scattering are generally larger than those produced in the radiation itself, particularly when 
the electron is low in energy. In spite of the possibly large angles between the electron and 
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the photon, the probability for this to occur is negligible and it can safely be concluded that 
bremsstrahlung and multiple scattering have no effect on the measured W boson mass. 

A last issue regarding radiation is the energy loss by ionization and by radiative processes 
where, for example, the electron radiates a photon that does not reach the calorimeter but 
produces an e + e~ pair that loses energy by ^P. These processes affect the W boson and 
Z boson mass and produce an offset in the energy scale, which was included in the energy 
scale determination. Small offsets produced in this way cancel to first order in the mean of 
the ratio of the W boson to Z boson masses, since the energy is lost to both the neutrino 
and the electron in each W boson event in which it occurs. In Z boson events only one 
electron loses the energy but the probability of such loss is twice as large. Using a GEANT 
simulation, a study of the effect of upstream material on the photon energy response was 
carried out. The photon response observed in the GEANT simulation was consistent with 
the response measured in situ, as described in Section [IV A| . Notably, the offset in response 
was found to be consistent with the in situ measurement. Combining all effects an overall 
systematic uncertainty of 20 MeV/c 2 was assigned to My/ due to radiative effects. 



K. Efficiency and Bias Uncertainties 

1. Trigger Efficiencies 

The effect of the uncertainty in the trigger efficiency has been studied by varying the 
nominal trigger efficiency distributions in the Monte Carlo within the range determined by 
the data. This resulted in an uncertainty on the W boson mass of 20, 20 and 60 MeV/c 2 
from the Mt, Pt an d Pr fits, respectively. In addition, the W boson mass was determined 
from a data sample that did not have the $t requirement imposed at the trigger level. The 
fitted mass from this sample was consistent with the nominal fit result within the statistical 
uncertainty, taking into account the large overlap between the two data samples. 
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2. Efficiency as a function of u» 

The transverse mass is relatively uncorrelated with the uncertainty in un, unlike the fits 
to the lepton transverse momentum spectra, which are very sensitive to this efficiency The 
nominal variation in the electron identification efficiency encompasses the band shown in 



Fig. 16. The results of large statistics Monte Carlo data samples generated with the nominal 



variations of the efficiency are given in Table |X11J|. Also listed are the results of the change 



in mass when fitting the data to templates generated with the different efficiencies. It is seen 
that the Monte Carlo studies and the data exhibit the same behavior. The corresponding 
electron identification uncertainty on the W boson mass is 20, 70 and 115 MeV/c 2 from the 
Mr, Pt and p v T fits, respectively. 



L. Error in the Fitting Procedure 

The W boson mass was obtained from an unbinned maximum likelihood fit in which the 
data were fit to transverse mass spectra which were generated for 21 different values of the 
W boson mass. The log-likelihood values for the different vector boson masses were fit to 
a parabola and the minimum was taken to be the fitted mass value. A decrease of half a 
unit in the log-likelihood was the one standard deviation statistical error. The likelihood 
distribution need not be Gaussian, depending on the range of the parameter fit, the intrinsic 
shape of the spectrum and the resolution function. The resulting log-likelihood curve was 
then non-quadratic. In addition, there will be fluctuations in the log-likelihood reflecting 
the Monte Carlo statistics. In order to determine the uncertainty, the fitting was redone 
with a cubic polynomial parameterization and the mass spacing was altered. This led to the 
assignment of 5 MeV/c 2 for the uncertainty due to the Mt fitting procedure. 
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M. Results of Systematic Errors 



The systematic errors on the W boson mass as obtained from the transverse mass, 
electron transverse momentum and neutrino transverse momentum are summarized in Ta- 
ble [XI V|. The measured mass results from this analysis are: 



M W (M T ) = 80.350 ± 0.140 ± 0.165 ± 0.160 GeV/c 2 ; 
M w (p e T ) = 80.300 ± 0.190 ± 0.180 ± 0.160 GeV/c 2 ; 
M w {f T ) = 80.045 ± 0.260 ± 0.305 ± 0.160 GeV/c 2 ; 



and from the transverse mass analysis, 

M w /M z = 0.88114 ± 0.00154 ± 0.00181 ± 0.00175. 

In each result, the first uncertainty is due to statistics, the second is due to systematic 
effects, and the third is due to the electron energy scale determination. 

VIII. CONSISTENCY CHECKS 

To verify the stability of the W boson mass result, consistency checks have been per- 
formed in which the W boson mass is determined from various modified data samples. These 
samples include those in which the fitting window was varied, additional selection criteria 
were applied, and a different electron clustering algorithm was used. Fully overlapping data 
samples were used to check the consistency of the results obtained from fits to the p 1 ^ and p v T 
spectra. Also, two-dimensional fits were done to check the consistency of parameters used 
in the Monte Carlo simulation. 

In general, the data sample was reduced or enlarged in these consistency studies. There 
was a large overlap between the original data sample and the samples used to verify the 
result. In order to quantify this verification, define the mass from the original data sample 
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M^? m , and that from the sample used in the verification M^ n . Then the estimator of the 
independent statistical error on the difference in the two results that were used is a(M^ m — 
M$ n ) = cr JEt. Here a is the statistical error on the original data sample, consisting of 
N\ + N2 events. The sample used for the consistency check contained N\ events. This is the 
error that is quoted for the difference in mass for the consistency checks. 



A. Additional Selection and Fitting Criteria 

To investigate the effect of multiple interactions, events were selected with low hit multi- 
plicity and a narrow time distribution in the small-angle scintillation counters (see Section 
II D). This yielded a sample in which approximately 77% were single interaction events. 
Also events with one and only one reconstructed event vertex and with only one track from 
the central detector in the electron road were selected. The latter cut removed mainly events 
with a random track from the underlying event. The change in fitted W boson mass from 
the transverse mass spectrum, with respect to the nominal mass value for each of these cross 
checks, is listed in Table |XV|. Note that the errors are statistical only. Any systematic error 



on the shifts is not included. 

To test the event modeling, the cut was tightened to 10 GeV/c and the result is 
listed in Table |XV| . When requiring p^ to be less than 10 GeV/c, there is an additional 



uncertainty due to the error on the hadronic energy scale factor and change in background 
contribution, which have not been included in the error estimate. 

Another important check of the event modeling is testing both the sensitivity and con- 
sistency of the result by tracking the change in mass during the process of applying different 
cuts. As an example, the first column in Table |XVI| lists the change in W boson mass from 



the nominal fit when u\\ in the data is required to be less than 10 GeV without modifying 
the templates. The change was rather dramatic for the mass from the p v T spectrum. The 
second column lists the change in mass when the templates are made consistent with the 
data. Even though the W boson mass is rather sensitive to the cut on u\\, the fitted masses 
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agreed well with the nominal values when data and Monte Carlo were treated consistently, 
indicating that both the scale and un efficiency were modeled correctly. 

To check for any systematic bias in detector response, event samples were selected with 
different fiducial requirements. For W boson events the rj range of the electron was restricted 
to electrons produced in the central region. For Z boson events the restriction was placed on 
only one of the two electrons, ensuring that a variation in detector response to the electron 



in W boson events was tracked in an identical manner in Z boson events. Table XVII lists 



the resultant change in the ratio of masses. The errors on the change were again statistical 



errors only. The ratio, tabulated in Table XVII, was with respect to the normalized ratio for 



the nominal 77 range. The ratio of masses did not change within errors. When the restriction 
was placed on both electrons in Z boson events, the ratio also did not change but had a 
large statistical uncertainty due to the significant loss of events. 

The variation in mass was also tracked when the nominal fitting range in transverse mass 



was varied. Figure |33| shows the change in W boson mass when varying the lower and upper 
edge of the fitting window for the fit to the distribution. Changing the fitting window 
led to a negligible systematic trend. 

B. Modified Electron Clustering 

Electron clusters were found by the reconstruction program using a nearest neighbor 
clustering algorithm The number of calorimeter towers included in the cluster was 

dynamic and depended on the environment of the electron. This algorithm thus introduced 
an uncertainty on the amount of underlying event energy included in the electron energy 
cluster, and therefore an uncertainty on how much energy was excluded from the underlying 
event for the calculation of the . In the discussion in Section |V C| the energy assignments 



and the modeling of the underlying event are described using a window algorithm for the 
reconstruction of the electron energy. The corrections necessary to translate these results 
to the cluster algorithm then had to be dealt with properly. These ambiguities can be com- 
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pletely circumvented if a fixed electron definition is used. To verify the internal consistency, 
the W boson mass was also determined using a fixed size electron cluster. 

The definition employed for the fixed size cluster is the "5x5 window algorithm". In 
this procedure, the electron energy was defined as the energy in the 25 towers in the region 
±0.2 in 7] and (p from the most energetic tower of the electron cluster as found by the nearest 
neighbor algorithm. Using the original W — > ev data sample, the electron energies were 
recalculated using the 5x5 window algorithm and the p^ was calculated with respect to the 
electron vertex, excluding the 5x5 window occupied by the electron. The region excluded 
from the underlying event for the calculation of the ]fi T was thus exactly known for each 
event. Subjecting these events to the standard event selection criteria yielded 7167 events, 
7131 of which were in the nominal data sample. The fitted W boson mass obtained from 
this data sample, using the window algorithm to define the electron, was 12 MeV/c 2 lower 
than when using the nearest neighbor algorithm. As noted above, a systematic uncertainty 
on the W boson mass of 20 MeV/c 2 has been attributed due to the difference in these two 
approaches and has been included in the underlying event uncertainty in Table (XTV. 



C. Fully Overlapping Data Samples 

The nominal fits to obtain the W boson mass were performed using events within a 
certain range either in transverse mass or in transverse momentum. These event samples 
did not fully overlap. Fully overlapping event samples are obtained when applying a fitting 
window in one variable and then utilizing the full unrestricted spectra in the other two 



variables, using all events in this window. Figure [H] shows the pf* and p? spectra with only 
the requirement that 60 < Mr < 90 GeV/c 2 . The change in W boson mass obtained from a 
fit to these spectra is +84 ±55 MeV/c 2 from the fit to the p e T spectrum and ±54 ±81 MeV/c 2 
from fitting the p v T spectrum. The errors on the shift in mass are the statistical errors due to 
the different number of events fit. Again, the results are consistent with the nominal results 
for the fits to the transverse momentum spectra. 
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D. Two-Dimensional Fits 



Two-dimensional fits were carried out to check the stability and correctness of parameters 
used in the Monte Carlo simulation. The first two-dimensional fit was performed on the W 
boson mass and the constant term in the electron energy resolution. Rather than expressing 
the likelihood in terms of the constant term, which resulted in a very asymmetric likelihood 
distribution, it was expressed in terms of the energy resolution at an electron pr of 40 GeV/c 



^40 = \/C2 + ^ (55) 

where S and C are the coefficients of the sampling and constant term, respectively. The 
sampling term was taken to be 0.13 and the constant term is varied. The error matrix for 
the fit in M w and TZ i0 is: 



( 0.0243 -0.0286^ 



y -0.0286 0.1933 j 



(56) 



with a correlation coefficient of p = —0.4155 . Figure [45] shows the contour in My/ and 
for a change of 0.5 units in the log-likelihood. The values on the axes are with respect to 
the central value of the fit. The fitted W boson mass was higher by 26 MeV/c 2 compared 
to the value obtained when the constant term in the energy resolution was fixed at 1.5%, 
in agreement with the nominal fit. The error on the mass from the two-dimensional fit was 
156 MeV/c 2 . For a fixed value of the constant term the error would be 0.156 a/1 — p' 2 = 
142 MeV/c 2 , consistent with the nominal fit result. The fitted value for the resolution 
(see Fig. f46]a) was 7Z i0 = (2.34 ± 0.440)% which, assuming a sampling term of S — 0.13, 
corresponded to a constant term of C = (1.1 jl°;f )%. This is again consistent with the result 
obtained from fitting the width of the Z boson resonance from which the constraint on the 
resolution is actually slightly tighter. The correlation between the W boson mass and 
was given by p = —147 MeV/c 2 per 1% change in 72-40, which was also consistent 

within errors with the result obtained in Section [VII C| and shown in Fig. |46|(b). 
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A two-dimensional fit was also performed in mass and hadronic energy scale factor k. 
The error matrix for this fit in My/ and K is: 



( 0.0250 0.0043 ^ 



^0.0043 0.003121 f 



(57) 



with a correlation coefficient of p = 0.457 . Figure |47] shows the one a contour in My/ and k. 
The values on the axes are with respect to the central value of the fit. The fitted W boson 
mass was lower by 7 MeV/c 2 compared to the value obtained when the pj? scale was fixed at 
0.83, in agreement with the nominal fit. The error on the mass from the two-dimensional fit 
is consistent with the nominal fit keeping the hadronic energy scale factor fixed. The fitted 
value for the hadronic energy scale factor was k = 0.834 ± 0.056, consistent with the result 
obtained from the Z boson data. The error is large because the W boson mass was not very 
sensitive to the hadronic energy scale. The correlation between the W boson mass and k is 
given by p = 12.8 MeV/c 2 per 1% change in scale factor. This is to be compared to 

the sensitivity of 12.1 MeV/c 2 per 1% change in scale factor obtained in Section [IV B| . 

In conclusion, the mass values obtained for different subsamples of the nominal data 
sample were all consistent within the quoted statistical uncertainty. Moreover, when leaving 
crucial parameters in the event modeling as free parameters in the fit, the W boson data 
preferred values for these parameters which were completely consistent with those obtained 
from external constraints, a strong indication of the stability of the result. 



IX. CONCLUSION 

A measurement of the W boson mass determined from the transverse mass distribution 
using electrons in the central region of the D0 detector from the 1992-1993 Fermilab Teva- 
tron running in 12.8 pb -1 has been described. The determination of My/ was based on a 
ratio of the measured W boson and Z boson masses, normalized to the world average Z 
boson mass as determined by the LEP experiments. This measurement yielded a W boson 
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mass value of 

M w = 80.350 ± 0.140 ± 0.165 ± 0.160 GeV/c 2 . 

and has an uncertainty comparable to that of other recent measurements in a single channel. 
The first uncertainty is due to statistics, the second is due to systematic effects, and the 
third is due to the electron energy scale determination. The 160 MeV/c 2 uncertainty due to 
the uncertainty on the absolute energy scale has a contribution of 150 MeV/c 2 due to the 
limited Z boson statistics. The measured ratio of the W boson and Z boson masses is 

M w /M z = 0.88114 ± 0.00154 ± 0.00181 ± 0.00175. 

Here, the first uncertainty is due to statistics, the second is due to systematic effects, and 
the third is due to the electron energy scale determination. 

Based on this measurement alone, the values for Ar and Ar res , as defined in Eqs. (0,f|), 
were determined to be: 

Ar = 0.0372 ± 0.0160 and 
Ar res = -0.0236 ± 0.0170. 



This measurement alone is thus sensitive to quantum corrections in the Standard Model at 
the 2.3cr level with evidence for bosonic radiative corrections with a significance of 1.4<r. 

An average W boson mass can be determined by combining the current result with recent 
previous measurements. The measurements are weighted with their uncorrelated uncertain- 
ties. The correlated uncertainty for the most recent measurements is that due to proton 
structure as parameterized in global parton distribution function fits. For each measurement, 
the uncertainty due to the common effect is removed to determine the uncorrelated error. 
Based on the UA2 [12|] and most recent CDF publication ITS] , the common uncertainty is 



taken to be 85 MeV/c 2 , the largest of the individual uncertainties due to the uncertainty on 
the structure of the proton. This procedure then yields a world average W boson mass of 
M w = 80.34 ±0.15 GeV/c 2 . 
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Figure ^ (top) presents a comparison of the world's direct determinations of My/ includ- 
ing this measurement and the overall pp world average. Also shown (band) is the Standard 
Model prediction using the LEP data as calculated by the LEP Electroweak Working Group 
fl9|| . Figure f|8] (bottom) shows the recently measured top quark mass Jl| from the D0 
collaboration versus the world average W mass. The top quark mass value used is 

m t (D0 ) = 172.0 ± 5.1 (stat) ± 5.5 (sys) Gev/c 2 (58) 

which is from the combined measurement of the dilepton and lepton plus jets channels. The 
Standard Model prediction for different values of the Higgs mass |55[] is also shown as the 
colored bands. 

Using the world average W boson mass, the derived values for the quantum corrections 
in the SM are 

Ar = 0.03834 ± 0.00885 and 
Ar rcs = -0.0224 ± 0.00944 . 

The direct measurement of the W boson mass at pp colliders indicates the existence of 
radiative corrections in the Standard Model at the ~ 4.3cr level and evidence of bosonic 
radiative corrections at the ~ 2.4<r level. 
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APPENDIX A: W AND Z BOSON PRODUCTION MODEL 

The theory and phenomenology of production of W and Z bosons can be divided into 
three regions of the pt of the vector boson. These regions are imprecisely ordered as follows: 

1. The high-j^ region in which perturbation theory is expected to be valid. This region 
is roughly 50 GeV/c and above. 

2. The low-pr region where perturbation theory is not helpful and soft gluons are freely 
emitted. There is a model for this process, and the validity of this theory is roughly 
below 15 GeV/c. By far, the bulk of the cross section for W boson and Z boson 
production is in this region. 

3. The intermediate region for which there is no theoretical description. Some analyses 
attempt to smoothly connect the two regions, beyond that which occurs naturally by 
simply adding the cross sections from region 1 to those of region 2. 

The Monte Carlo generation of the vector bosons relied on the resummation formalism 
of Collins, Soper, and Sterman (CSS) P9| which treats the emission of soft gluons in region 
2 by summing all contributions in impact parameter space. There are few free parameters 
in this model and it is shown below that it satisfactorily matches the D0 data. The triple 
differential cross section for production of a W boson can be written 



Here, Wij(b*, Q, x A , %b) includes the convolution of parton densities for partons i,j and the 
splitting functions, the Cabibbo-Kobayashi-Maskawa elements, and the electroweak param- 
eters. The quantity Q is the invariant mass of the annihilating partons, while xa,b is the 
Bjorken x variable representing the fraction of the colliding hadron's momenta carried by 




^W ij (b*,Q,x A ,x B )e- s ^^F^ 



(b,Q,Q ,x A ,x B ) 



+Y(p T ,Q,x A ,x B )- 



(Al) 
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the annihilating partons. Qo is taken to be the lowest scale where perturbation theory is 
presumed to be sensible. The quantity o"o is for normalization. The Sudakov form factor 
S(b, Q) is fixed by the order in a s and is an integral over a running scale. The combination 
of these terms describe region 2. The quantity Y(pt,Q,xa,%b) contains terms which are 
less singular than Pt~ 2 and is the term which dominates in the perturbative regime, region 
1. 

The complication inherent in this formalism is the Fourier transformation of the cross 
section, which involves an integral over all values of the impact parameter b. This is dealt 
with by regulating b to behave well near the origin, forcing it to tend to a constant as 
b ^ 0. In the CSS formalism, this amounts to a replacement of b — > b* = -, b . The 
price for making this modification is the obligation to add a term to "replace" the missing 
contribution to the integral from this b — > b* substitution. This extra factor is the so-called 
non-perturbative function, represented in Eq. as F[- p . Theoretical arguments fix the 
form of F^ p , up to phenomenological parameters. 

There have been two efforts to determine the non-perturbative function. One such recent 
fit is by Ladinsky and Yuan El (LY) who parameterized the non-perturbative function as 



(A2) 



F^ p (b,Q,Q ,x A ,x B ) = 

Q 



cxp 



b 2 gi - g 2 b 2 \n I ) - 5-15-36 ln(100xAX B ) 



The g parameters are not specified by theory, but are measurable. A much earlier effort 
by Davies and Stirling EJ (DS) used an identical parameterization, but essentially with 
53 = 0. Recently, Arnold and Kaufman |^0| (AK) employed the CSS formalism including 
the DS fits, a NLO calculation for the Y term [^T| (region 1), and a strategy of dealing with 
region 3. A computer program has been available for the AK approach. Likewise, the LY 
calculation was done with an independent computer program which is identical in its coding 
of the CSS theory, but utilized a simple 0(as) calculation for the Y term. LY made no 
attempt to match the two regions. Both computer codes have been used in this analysis. 
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This description of W boson and Z boson production is taken as the ansatz for the 
Monte Carlo production model. The more recent LY fits to modern Drell-Yan and collider 
Z boson data constrain the g parameters and have been used here as representative of the 
best available information. In this sense the g parameters function operationally like the 
parton distribution functions. The LY fits result in g x = 0.1 GeV 2 , g 2 = 0.58tg;2o 
GeV 2 , and g% = —1.5 ± 0.10 GeV -1 . These central values have been used as the nominal 
production model for W and Z bosons, with the major sensitivity to c/2- 



APPENDIX B: BREMSSTRAHLUNG 



The fraction of decays which involves radiation depends on the minimum photon energy, 



E™ m , which was taken to be 50 MeV. Figure ^ shows this fraction as function of E™ m for 



7 



(a) W boson and (b) Z boson decays. For Z boson decays the fraction of radiative decays 
is about a factor of two higher than for W boson decays, as expected. For the default E™ m , 
31% of the W boson decays and 66% of the Z boson decays were radiative. Only order 
corrections have been included and so processes in which two or more photons are radiated 
were not generated. 

For radiative W boson decays, W — *> ewy, it is important to determine the minimum 
spatial separation between the photon and electron that would result in the photon energy 
not being included with that of the electron by the reconstruction program. For events with 
R = y/ Ai] 2 + A(/? 2 above approximately 0.2 the photon energy may not be added to that 
of the electron. Instead, it was reconstructed as part of the W boson recoil. The neutrino 
energy was unchanged, but the electron energy is too low. The W boson and Z boson masses 
were then too low in a manner which does not cancel in the ratio. Since this effect is large, 
it is important to evaluate the effect when the photon is produced by bremsstrahlung in the 
central detector. 

For the photon to have an effect on the measured W boson mass, it must be separated 
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from the electron in (77, ip) space by at least R ei = 0.2, that is, 

A v 2 + Ap 2 > R 2 er (Bl) 
With An = -=^; = cosh ri A$, this can be written as 

' sin # ' ' 

2 

+ A^ 2 > R 2 (B2) 



y sin $ 

Switching to coordinates measured with respect to the electron 

At? = a; cos a (B3) 
= "Tf . (B4) 

where w is the angle between the electron and the photon and a the azimuthal angle of the 
photon with respect to the electron, one can write: 

(to cos a) 2 + (to sin a) 2 > ( — P—\ (B5) 

y cosh 77 J 

or 

u > (B6) 
cosh?^ 

The angle between the electron and photon must be greater than 0.2 rad for rj = and 
greater than 0.13 rad for 77 = 1. In units of where m e is the electron mass and E the 
electron energy, this corresponds to — x uj > 13, 000 for an electron energy E of 50 GeV 
and 7] = 1. 

Figure 0(a) shows the probability ^ for radiating a photon at an angle u for the case 
y=0.1, where y is the fraction of the electron's energy carried by the photon fffi] . The angle 
is expressed in units of For all calculations in this analysis, Z=13 (aluminum) has 
been assumed and the energy of the electron has been fixed to E = 50 GeV. The probability 



decreases by four orders of magnitude at ^-scaled angles of 50. Figure |50|(b) shows the 
relative probability for radiating a photon at an angle u and its dependence on y. Although 
the probability for radiating a photon is larger at small y, after normalization, there is little 



y dependence of the angle at which the photon is radiated. Since scaled angles of 13,000 
or more are needed, the photon never separates from the electron by radiation alone and 
therefore bremsstrahlung has no effect on the W boson mass. 

The electron and photon can also separate if the electron scatters through a large angle. 
The probability that an electron radiates a photon of momentum between k and k + dk in 
dx in a medium with radiation length Xo is [|53"l,[48j 

, \ „ , dx dk ( 4 4 9 \ , 
P(E,k)dkdx = —— -y + y 2 y (B7) 

The quantity y is the fraction of the electron's energy carried by the photon, y = J|, and E 
and k are the electron and photon energies. Integrating from k = k m i n to E one finds 

dx 4 

P(E,k > k min )dx = — [--(ln?/ min + 1 - y min ) 

A 6 

+ -vLn)]dx. (B8) 

For y min close to 1 this gives 

fj nr j 

P(E, k > k min )dx = — (C + -C 2 )dx (B9) 

X 6 

with ( = 1 — y m i n . For example, the probability that a 50 GeV electron radiates a 49 GeV 
photon in 0.15 Xo is 

0.15 x 0.02 « 3 • 10~ 3 . (B10) 



The 1 GeV electron can then scatter through a large angle JT8 

w rms = V2 x 13 ' 6 MeV x y/oJs = 7.4 mra d. (Bll) 
1 GeV 

The angles resulting from multiple scattering are generally larger than those produced in the 
radiation itself, particularly when the electron is low in energy (u Tms ~ 4). Nevertheless, it is 
still difficult to separate the photon and electron sufficiently The 7.4 mrad angle calculated 
above translates, in units of into — x 7.4 x 10~ 3 ~ 800, still small compared to 13,000. 
The falloff in scattering is rapid (Gaussian). If one considers a 50 GeV electron radiating 
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99.8% of its energy, the probability becomes small, 3 x 10~ 4 . The resulting 100 MeV electron, 
however, can now multiple scatter through a large angle, 80 mrad or 8000 in units of 

This situation is compared with radiative W boson decays in Fig. |5"I|(a) which shows 
the distribution in uj, in units of ^ , for radiative W boson events with the electron in the 
central calorimeter with > 25 GeV/c. The distribution has a very long tail extending to 
values of 50,000 for uj. At small angles of ^- x uj = 10, 000 the cross section is down by a 
factor of approximately 200. Nevertheless, 21% of the events have angles greater than 5000. 
Figure |5l](b) shows the event distribution in uj for events in which the photon and electron 
reconstruct as separate entities with the photon retaining its identity. As was estimated 
above, the threshold at x uj is approximately 15,000. 

APPENDIX C: MEAN NUMBER OF INTERACTIONS 

The library of minimum bias events was stored in bins of luminosity according to the 
following rule. Given a W boson event, recorded at a luminosity £ with corresponding 
average number of interactions per crossing (n), the minimum bias event, mimicking the 
underlying event, was taken at a scaled value of the instantaneous luminosity, £'. £ was 
chosen so that the mean multiple interaction rate in Monte Carlo generated W boson events 
is the same as in the W boson data sample. 

The probability of getting a W boson trigger in a crossing in which there are n interactions 
is given by 

P(W,n) = nP(n)— . (CI) 

0"inel 

Here Pin) is the Poisson probability of n interactions in the crossing, ■ £J£ - the probability 

CT inel 

that the inelastic interaction is one in which a W boson is produced. The factor n represents 
the number of ways one can choose the W boson interaction from the n interactions in the 
crossing. Note that the probability of getting a W boson in a crossing is then 

P(W) = ]T nP(n)^ = (n) (C2) 

n 0~inel 0"incl 
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which is the expected rate when (n) is written in terms of the luminosity and the inelastic 
cross section, Oi ne i. The probability distribution of getting n interactions in a crossing in 
which there is a W boson is 



P{n\W) = nP(n), (C3) 

and has a mean value of (n) + 1. This shows that the mean number of interactions in a 
crossing in which there is a W is (n) + 1. 

For the minimum bias trigger, the average number of interactions per crossing (n m i n ) 
given that there is at least one, is 

_ SStpn'PCn') _ (n>) 
\n mm) Y%=iP(ri) 1-e-W 1 1 

The minimum bias events are chosen at a luminosity £' such that the mean (n TO j„), as 
given by equation |C4], is equal to (n) + 1, where (n) is the mean number of interactions at 



luminosity L at which the W event was recorded. This guaranteed that the mean number 
of interactions was correct. The distributions in the number of interactions per crossing 
are somewhat different, though. The minimum bias distribution is a Poisson distribution, 
cut off at n = 1, while the number of interactions in W events is a Poisson distribution, 
beginning at n = 1. The impact of this difference in this analysis is negligible. 
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FIG. 1. Kinematic quantities for W events. 
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FIG. 2. (a) Transverse mass distribution of W events and (b) dielectron invariant mass distri- 
bution from Z events. Both distributions are shown before the application of fitting window cuts 
or energy scaling (see Section IV). 
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FIG. 3. Observed invariant mass from Z — ► ee decays, m™ eas , versus 
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FIG. 4. Constraints on a and 5 from a) Z — > ee decays (solid contour), b) J/ip —> ee decays 
(dotted lines), and c) ir° —>■ 77 — ► 4e decays, (dashed lines). 
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FIG. 5. The invariant mass from 7r° — > 77 — > e + e e + e decay events (points). Also shown is 
the background contribution (open circles). 
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FIG. 6. Background subtracted invariant mass from the 7r° event sample (points) compared to 
the Monte Carlo simulation (line). 
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FIG. 7. Mass distribution of the observed J /ip — > ee decays (histogram) is shown above the 
background (points). The line is a fit to the signal plus background. 
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FIG. 8. Expanded view of Fig. [I] showing a versus 5 with the combined best fit (shaded region). 
The expanded lobe (dotted contour) to lower values of 5 is due to uncertainties in the low energy 
non-linear response of the calorimeter. The contributions are from: Z — > ee decays (solid contour), 
J/ip — > ee decays (dashed-dotted lines), and ir° — > 77 — > e + e~e + e~~ decays (dashed lines). 
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FIG. 9. Definition of the 7?-£ coordinate system for Z events. 
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FIG. 10. For Z — ► ee events (points) (a) the average r\ imbalance versus pj? • f) is shown along 
with the line obtained from a linear least-squared fit to the data and (b) the rj imbalance with a 
hadronic energy scale factor k = 0.83 applied is shown with a Gaussian fit (curve). 
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FIG. 11. For Z — > ee events (points) with the same event topology as W events, the average 
value of pj? c • V is shown versus p^ e ■ fj. The line shown is obtained from a linear least-squares fit 
to the data. 
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FIG. 12. Distribution of the average + + $t)\ 2 versus |pf 6 | 2 for Z — ► ee events. The 
line shown is obtained from a linear least-squares fit to the data. 
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FIG. 13. (a) The modeled resolution of the z position of the center of gravity of CDC tracks 
and (b) the distribution in the difference of the intersections of the z axis of the two electron tracks 
from Z decays (points) compared with the distribution from the Monte Carlo simulation. 
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FIG. 14. Average value of the isolation versus for electrons from W ^ ev decays. 
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FIG. 15. Distribution of the isolation value, fi so , for electrons from W — > ev decays for different 
ranges (points). The curves are fits to the data. 
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FIG. 16. Electron identification efficiency as a function of u» (open crosses) from data. The 
central curve is a fit to the data. The outer curves show the allowed ranges for determining the 
systematic errors. 
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FIG. 17. Transverse mass spectrum of the multi-jet background obtained from the data (open 
crosses). The solid line is a fourth-order polynomial fit. 
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FIG. 18. Spectra in (a) My, (b) an d (c) $r for the Z boson background in the W boson 
sample. The lines are fits to the data. 
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FIG. 19. Normalized distributions in transverse mass of the dominant background contributions 
to the W boson event sample. 
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FIG. 20. Distribution in instantaneous luminosity of the W events used in the W boson mass 
measurement. 
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FIG. 21. (a) ^ and (b) T,Et distributions of the minimum bias events used to model 
W — > eu and Z — ► ee underlying event. 
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FIG. 22. Sensitivity of the width of the r\ imbalance distribution to the number of minimum 

bias events used to simulate the underlying event in the Monte Carlo simulation (points). The line 
is the result of a linear least-squares fit. The bands (dotted-dashed) correspond to the nominal 
and ± la measurements of the width in Z boson events. 
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FIG. 23. Angular distribution of electrons from W — > eu decays (points) compared to the 
simulation (histogram). The asymmetry is due to the fact that the luminous region was not 
located at z = cm in the D0 detector, but was rather centered at z = —7.98 cm. 
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FIG. 24. Distribution of from W — ► ev decays (points) compared to the simulation (his 
gram) . 
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FIG. 25. Distribution of the angle between the recoil jet and the electron in the transverse 
plane from W — ► eu decays (points) compared to the simulation (histogram). 
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FIG. 26. Distribution of the mean un versus p?p from — > ev decays (points) compared to the 
simulation (*). 
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FIG. 27. Distribution of the mean uu versus pj- from W — > ev decays (points) compared to the 
simulation (*). 
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FIG. 28. Distribution of the mean uy versus Mt from W — * ev decays (points) compared to 
the simulation (*). 
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FIG. 29. Distribution of the mean uii versus from W — > ev decays (points) compared to 
the simulation (*). 
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FIG. 30. Comparison of the u\\ distribution from W — > eu events (points) and the Monte Carlo 
simulation (histogram) . 



119 



> 




u ± (GeV) 



FIG. 31. Comparison of the u± distribution from W — ► ev events (points) and the Monte Carlo 
simulation (histogram) . 
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FIG. 32. Comparison of the electron energy distribution from Z — > ee events (points) and the 
simulation (histogram). 
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FIG. 33. Comparison between the Z boson transverse momentum distribution as measured 
from the recoil system in Z — > ee events (points) and the simulation (histogram). 
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FIG. 34. Distribution of (a) the fitted masses and (b) the fit uncertainties from fits to the 
transverse mass distributions for an ensemble of 125 Monte Carlo generated data samples of 8000 
W — ► ev decays. 
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FIG. 35. (a) The central dielectron invariant mass distribution for Z events (points) and the 
best fit of the simulation (histogram), (b) the corresponding relative log- likelihood distribution and 
(c) signed \fx* distribution. 
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FIG. 36. (a) The transverse mass distribution for W events (points) and the best fit of the 
simulation (histogram), (b) the corresponding relative log- likelihood distribution and (c) signed 
\fx? distribution. 
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FIG. 37. (a) The electron transverse momentum distribution for W events (points) and the 
best fit of the simulation (histogram), (b) the corresponding relative log- likelihood distribution 
and (c) signed v5? distribution. 
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FIG. 38. (a) The neutrino transverse momentum distribution for W events (points) and the 
best fit of the simulation (histogram), (b) the corresponding relative log- likelihood distribution and 
(c) signed distribution. 
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FIG. 39. Average transverse energy flow per electron cluster tower as a function of r\ measured 
from minimum bias events. 
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FIG. 40. Ratio of predicted differential cross section in and the nominal cross section for 
new parameterizations of the CTEQ3M parton distribution function. 
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FIG. 41. (a) Ratio of predicted differential cross section in and the nominal cross section 
versus p^ when the parameter <?2 is varied by multiple standard deviations from its nominal value 
in the Ladinsky-Yuan prediction and (b) the distribution in x 2 f° r a comparison between data and 
Monte Carlo of the p^, spectrum versus the variation of g2 in units of its standard deviation. 
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FIG. 42. The measured multi-jet background distribution versus My from the data (open 
crosses). The allowed variations in the shape of the transverse mass spectrum (dotted lines) are 
shown. The solid line indicates the nominal background distribution. 
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FIG. 43. Change in fitted W boson mass when varying the (a) upper and (b) lower edge of 
the fitting window from the fit to the transverse mass spectrum (points). The horizontal bands 
indicate the la statistical error on the nominal fit. 
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FIG. 44. The (a) electron and (b) neutrino transverse momentum distribution for the events 
in the transverse mass window 60 < Mt < 90 GeV/c 2 (points). The histograms are the best fits 
of the simulation. 
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FIG. 45. The la contour in the change in M\y and the electron energy resolution at pj, = 
40 GeV/c from fits of the simulation, in which the constant term is allowed to vary, to W boson 
events. 
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FIG. 46. From fits of the simulation, in which the constant term is allowed to vary, to W events: 
(a) The relative likelihood versus the electron energy resolution at pf. = 40 GeV/c (b) The change 
in the fitted W mass versus the electron energy resolution at pj, = 40 GeV/c. 
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FIG. 47. The la contour in the change in My/ and scale from fits of the simulation, in 
which the scale factor is allowed to vary, to W boson events. 
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FIG. 48. The upper half of the figure shows the D0 determination of Mw along with recent 
results from other hadron collider experiments and the pp world average (see the text for a discussion 
of the world average calculation). The band is the Standard Model prediction from the combined 
LEP results. The lower half of the figure shows the D0 determination of the mass of the top quark 
versus the world average determination of Mw (*)■ The contour shows the allowed range in each 
value. The Standard Model prediction (see the text) for various assumptions of the Higgs boson 
mass is indicated by the bands. 
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FIG. 49. Fraction of radiative (a) W decays and (b) Z decays as function of the minimum 
photon energy. 
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FIG. 50. (a) Probability for an electron to radiate a photon in aluminum (Z=13) as function 
of the angle u between the electron and the photon in units of ^f, where m e is the electron mass 
and E its energy, (b) Relative probability for radiating a photon for different values of y = -g, 
with k the photon energy. 



139 



Radiative decays central electrons 



Central Electrons 

p T (e) > 25 GeV/c 
E min ( T ) = 50MeV 




10000 20000 30000 40000 50000 

co (m e /EJ 



Radiative decays central electrons 



o 
o 
o 



0) 

c 
a 

&io 2 



10 



Photon Identity Preserved 



Central Electrons 

p T (e) > 25 GeV/c 
E mi „(Y) = 50 MeV 




10000 20000 30000 40000 50000 

to (m./EJ 



FIG. 51. Distribution in u, in units of for radiative W boson events in the central calorime- 



ter and (b) for events where the the photon retains its identity. 
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TABLES 



Experiment 


channel 


Mw (GeV/c 2 ) 

-i rJ. yy \ VJ v i / v / 


UAlflQS.^ H 

yj -ii- ±. \ ±. (j : ) i \y~'\\ 




81 ± 5 


UA2flQ83 N ) Hi 




80+^° 


TJA1('1Q86 N ) Hi 




83 5 + H ±27 


UA2fl987 N ) H 

\_/ .a 1. — ( \ J. kj i / l t -'l 




80 2 ± 6 ± 5 ± 1 3 

v. J V / • 1 V J • \J 1 V J ■ '.J 1 _L • (J 


TjAifiQsq'i [Full 




82 7±lf)±27 




ez/ 


an /1Q 4- n /iq 4- n O/i 

oU. 4y ± U.4o ± U.Z4 


UA2fl992) [O] 


ev 


80 35 ± 33 ± 17 


CDF(1989) pfl 


ev 


80.0 ± 3.3 ± 2.4 


CDF(1990) El 


cv 


79.91 ± 0.35 ± 0.24 ± 0.19 


our^iyyoj ||10|| 


cv 


Sfl AQfl 4-014^4-017^ 


u-rii ^lyo^j ||io|| 


[IV 


81 + 6 

Ol_y 


UA1(1989) p[| 


[IV 


81.8 + f-Q ± 2.6 


CDF(1990) 0I 


[IV 


79.90 ± 0.53 ± 0.32 ± 0.08 


CDF(1995) || 


[IV 


80.310 ±0.205 ±0.130 


UA1(1989) @ 


TV 


89 ± 3 ± 6 


CDF(1990) 


cv + [iv 


79.91 ± 0.39 


CDF(1995) || 


ev + [iv 


80.410 ±0.180 



TABLE I. Previously published hadron collider measurements of Mw. In each case the first 
uncertainty listed is statistical, the second is systematic, and the third is due to energy scale. 
For the latest CDF values, the energy scales have been incorporated into the total systematic 
uncertainty. 
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W — > ev candidates Z — > ee candidates 



LI trigger requirements 1 EM tower with £r > 10 GeV 2 EM towers with Er > 7 GeV 
L2 filter requirements 1 EM cluster with E T > 20 GeV 2 EM clusters with E T > 10 GeV 

fiso < 0.15 / iso < 0.15 

$ T > 20 GeV 

TABLE II. LI and L2 trigger requirements for W and Z event data samples. Here, fi so is 
defined by Eq. 



W boson event sample 



ECN 
1838 


CC 
7234 




ECS 
1681 






Z boson event sample 
ECN-ECN ECN-CC CC-CC 


CC-ECS 


ECS-ECS 


Mass Measurement 
Resolution Studies 


48 147 366 
46 143 344 


134 
130 


39 
35 



TABLE III. Event samples for W and Z bosons. Here, "N" and "S" refer to the end calorime- 
ters on the north and the south. 
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Descriptor 



Nominal value 



EM energy resolution, sampling (CC) 
EM energy resolution, constant (CC) 
EM energy resolution, noise (CC) 
HAD energy resolution, sampling (CC) 
HAD energy resolution, constant (CC) 
HAD energy resolution, noise (CC) 
HAD energy scale 
Electron Underlying Event 
W Width 
Z Width 

# minimum bias events 



minimum E. 



7 



AR(e-f) 

Calorimeter position resolution 



CDC z cog resolution 



if resolution 



S = 13.0% 
C = 1.5% 
N = 0.4 GeV 
S = 80.0% 
C = 4.0% 
N = 1.5 GeV 
k = 0.83 
Ef} E = 205 MeV 
IV = 2.1 GeV 
T z = 2.5 GeV 
1.0 

E min = 50 Me y 

AR(e-f) = 0.3 
a(z) w 0.7 cm 



rms z cog = 2.0 cm 



a(<p) = 0.005 rad 



TABLE IV. Parameters used in the fast Monte Carlos. 
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W production 



Z° production 



pdf 


(3 x 100 
(GeV- 1 ) 


v — v and 

(%) 


v — s s — s 

(%) 


(3 x 100 v - i 
(GeV- 1 ) 


' and v — s 

(%) 


s — s 
(%) 


MRS E' 


0.980 


82.7 


17.3 


0.869 


84.7 


15.3 


MRS B 


1.054 


82.7 


17.3 


0.897 


85.0 


15.0 


HMRS B 


1.048 


75.5 


24.5 


0.932 


77.7 


22.3 


KMRS B0 


1.022 


79.2 


20.8 


0.908 


81.4 


18.6 


MRS DO ' 


1.220 


78.9 


21.1 


1.077 


80.9 


19.1 


MRS D'- 


1.277 


79.9 


20.1 


1.097 


81.7 


18.3 


MRS H 


1.264 


79.0 


21.0 


1.104 


81.0 


19.0 


MRS A 


1.282 


79.6 


20.4 


1.101 


81.0 


19.0 


MRS G 


1.297 


80.3 


19.7 


1.107 


81.6 


18.4 


MT Bl 


1.076 


83.1 


16.9 


0.925 


85.4 


14.6 


CTEQ 1M 


1.204 


79.6 


20.4 


1.038 


81.3 


18.7 


CTEQ IMS 


1.206 


79.9 


20.1 


1.030 


81.6 


18.4 


CTEQ 2M 


1.274 


79.4 


20.6 


1.078 


81.0 


19.0 


CTEQ 2MS 


1.231 


79.7 


20.3 


1.043 


81.2 


18.8 


CTEQ 2MF 


1.225 


78.7 


21.3 


1.054 


80.2 


19.8 


CTEQ 2ML 


1.310 


79.7 


20.3 


1.113 


81.4 


18.6 


CTEQ 3M 


1.224 


79.7 


20.3 


1.051 


81.1 


18.9 


GRV H0 


1.237 


82.0 


18.0 


1.095 


80.5 


19.5 


TABLE V. Parton luminosity slope, valence- valence (v — v), valence-sea 
(s — s) contributions to the W and Z boson production cross section at \fs= 


(v — s) and 
1.8 TeV. 


sea-sea 
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rj response 


AM W 


AM W 


AM W 


AM Z 




M T fit (MeV/c 2 ) 


p e T fit (MeV/c 2 ) 


p u T fit (MeV/c 2 ) 


m ee fit (MeV/c 2 ) 


module A 


-6 ± 16 


-7 ± 22 


-49 ± 30 


-2 ± 6 


module B 


+5 ± 16 


-15 ± 22 


-26 ± 30 


-8 ± 6 



TABLE VI. Change in W and Z boson masses in MeV/c 2 if a non-uniform calorimeter 77 
response is assumed, bracketed by the variations observed for two EM modules exposed in a test 
beam. 



fitted spectrum Monte Carlo Sensitivity Data 

9Mvi 
~5U 



AM W (MeV/c 2 ) (^v/£!) AM W (MeV/c 2 ) 



M T tu ± 17 -112 ± 19 +^ 



If ± 22 -54 ± 14 

p v T t%t ± 30 -56 ± 19 tf 

TABLE VII. Uncertainty on the II 7 boson mass in MeV/c 2 due to a change in the constant 
term of the electromagnetic energy resolution by 0.5%. The upper numbers correspond to the 
lower constant term. 
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fitted spectrum 


Monte Carlo 
AM W (MeV/c 2 ) 




Sensitivity 

dM w /MeV/c 2 N. 
dn V 0.01 > 




Data 
AM W (MeV/c 2 ) 


M T 
Pt 
Pt 


+55 i -i y 

+38 _i_ 90 
-29 31 zo 

;gf ±30 




+ 12.1 ± 1.3 

+6.7 ± 1.7 
-30.3 ± 2.5 




+42 
-80 

+4 
-38 

-125 
+100 


TABLE VIII. Uncertainty on the W boson 


mass due to the change 


w 
m prf, 


scale by 0.04. The 


upper numbers are 


the change in mass when the 


p¥ 


scale factor increases 


and the hadronic response 



is closer to the electromagnetic response. 



fitted spectrum Monte Carlo Sensitivity Data 

dM w 
d#min.bias 



AM W (MeV/c 2 ) (M£^) AM W (MeV/c 2 ) 



M T +m±17 -117 + 5 

~\l ± 23 -20.0 ± 7.0 ~f 



-245 _i_ qh oc« J- ~\ A -535 

+554 



± 30 -286 ± 14 



TABLE IX. Uncertainty on the W boson mass due to a change by 0.1 in the number of 
minimum bias events underlying the W event. The upper numbers are the change in mass for a 
higher average number of minimum bias events. 
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fitted spectrum 


Monte Carlo 


Sensitivity 


Data 




AM W (MeV/c 2 ) 


dM w /MeV/c 2 ^ 

(JO v 1U/0 ' 


AM W (MeV/c 2 ) 


lVl r p 


-74 117 
+52 ± 17 


— Ol. ± D.U 


-52 
+31 


Pt 


±| ± 23 


-2.5 ± 7.8 


—4 
-26 


P T 


i| ± 30 


-38.3 ± 11.0 


-87 
+32 


TABLE X. 


Uncertainty on the W boson mass 


due to the change in 


the sampling term of 


the hadronic energy resolution by 0.2 . The upper 


numbers are the chan^ 


;e in mass for a larger 



resolution. 



PDF (CTEQ3M) g 2 - 2a g 2 g 2 + 2a g 2 + 4cr fit 

AM W (MeV/c 2 ) AM W (MeV/c 2 ) AM W (MeV/c 2 ) AM W (MeV/c 2 ) 
CDF Asym. high +32 +14 +50 +11 M T 

CDF Asym. nominal -14 -37 —30 Mt 

CDF Asym. low -55 -67 -69 -65 M T 

CDF Asym. high +125 +51 +36 

CDF Asym. nominal +45 -93 

CDF Asym. low -48 -127 -169 

CDF Asym. high +64 +80 +77 

CDF Asym. nominal +40 -43 

CDF Asym. low -64 -69 -141 

TABLE XI. Shift in the W boson mass in MeV/c 2 when using different parametrizations of 
the parton distribution functions and pYp spectrum. There is a statistical uncertainty of 17, 24 and 
31 MeV/c 2 on each value for the Mp, Pj> and p v T fit, respectively. 



-60 
-137 
-197 



Pt 



Pt 



Pt 



-17 
-78 



Pt 



Pt 



-121 p T 
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PDF 


AM W 
M T fit (MeV/c 2 ) 


AM W 
p e T fit (MeV/c 2 ) 


AM W 

p u T fit (MeV/c 2 ) 


MRSA 


— 


— 


— 


mrsbM 


-90 ± 19 


-196 ± 24 


-86 ± 34 


mrseW 


-136 ± 19 


-168 ± 24 


-198 ± 34 


hmrsbW 


-157 ± 19 


-280 ± 24 


-204 ± 34 


KMRSB0M 


-175 ± 19 


-238 ± 24 


-244 ± 34 


MRSDO' 


-74 ± 19 


-109 ± 24 


-26 ± 34 


MRSD'- 


-31 ± 19 


-9 ± 24 


+8 ± 34 


MRSH 


-30 ± 19 


-47 ± 24 


-70 ± 34 


MTBlW 


-135 ± 19 


-260 ± 24 


-144 ± 34 


CTEQIMSM 


-29 ± 19 


-109 ± 24 


-1 ± 34 


CTEQ2M 


+20 ± 19 


+1 ± 24 


+53 ± 34 


CTEQ2MS 


± 19 


-26 ± 24 


+62 ± 34 


CTEQ2MF 


-59 ± 19 


-112 ± 24 


-84 ± 34 


CTEQ2ML 


+29 ± 19 


+19 ± 24 


+57 ± 34 


CTEQ3M 


-31 ± 19 


-75 ± 24 


-102 ± 34 


GRVH0 


-47 ± 19 


-88 ± 24 


-50 ± 34 


TABLE XII. Change 


in the W and Z boson 


masses in MeV/c 2 with varying 


; parametrizations 



of the structure of the proton for transverse momentum spectra. Amounts quoted are relative to 
the MRSA fit. The asterisk indicates those parton distribution functions considered obsolete for 
this analysis. 
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Monte Carlo 


Data 


fitted spectrum 


AM W (MeV/c 2 ) 


AM W (MeV/c 2 ) 


M T 


tf ± 17 


+2 
-13 



e -46 , 2 o -63 

"T +52 ^ zo +41 

n v +124 , on +136 

PT -143 ^ JU -95 

TABLE XIII. Uncertainty on the W boson mass due to uncertainty on the electron identifica- 
tion efficiency as a function of the quantity uy. The upper numbers are the change in mass when 
the overall efficiency decreases. 
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Source 


Variation Used 


a(M w ) 


a(M w ) 


a(M w ) 






M T Fit 


p e T Fit 


p v T fit 






(MeV/c 2 


) (MeV/c 2 ) 


(MeV/c 2 ) 


Statistical 




140 


190 


260 


Energy Scale 




160 


160 


160 


Other Systematic Errors 




165 


180 


305 


EM energy resolution 


c = (i.5l?:i) 


70 


35 


35 


CDC z scaled 


a = (0.988 ± 0.002) 


50 


55 


55 


Hadronic energy resolution 


S had = 0.8 ± 0.2 


65 


5 


80 


Underlying event*-*) 


^Tower = ( 16 .g ± 1-5 ) MeV 


35 


35 


35 


T w 


T w = (2.1 ±0.1) GeV 


20 


20 


20 


Hadronic energy scale 


a had = (0.83 ± 0.04) 


50 


30 


120 


Number of minimum bias events (1.0 ± 0.06) 


60 


10 


150 


QCD background 


(1.6 ±0.8)% 


30 


35 


35 


Z — ► ee background 


(0.43 ± 0.05) % 


15 


20 


20 


Electron ID efficiency 


parametrization 


20 


70 


115 


Radiative decays 


Tpmin t> ,.2 


20 


40 


40 


p T (W), pdf 


Pt(W) variation 


65 


130 


130 


Trigger efficiency 


efficiency spread 


20 


20 


60 


Non-uniformity in r] 


test beam 


10 


10 


25 


Fitting error 




5 


10 


10 


Total 




275 


315 


435 


TABLE XIV. Summary 


of systematic errors on the W boson mass 


from the three 


mass fits. 



Those errors that are strongly correlated with the measured Z boson mass are indicated by an 
asterisk. 
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Data subsample 



AM W (MeV/c 2 ) 



One track in electron road^*) -2 ± 54 

One reconstructed event vertex^ -76 ± 76 

Single interaction events W -107 ± 95 

< 10 GeV/c -166 ± 90 

TABLE XV. Change in W mass from nominal for different subsamples of the data. Those 
subsamples for which the Monte Carlo templates were not modified are indicated by an asterisk *. 
Errors are statistical only. 



fitted spectrum AM W (MeV/c 2 ) AM W (MeV/c 2 ) 

Data, no u\\ cut Data, ity < 10 GeV 

MC, U || < 10 GeV MC, < 10 GeV 



M T 



Pt 



P u t 



+78 
-280 
+810 



-16 

+40 
-45 



TABLE XVI. Change in W mass from nominal when applying a cut on u\\ of 10 GeV. 
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r\ range R = 



M z 



\r] e \ < 1.0 1.0003 ± 0.0005 

\rj e \ < 0.8 1.0010 ± 0.0012 

\rj e \ < 0.6 1.0011 ± 0.0019 

TABLE XVII. The ratio of the W and Z boson masses when restricting the rj range of the 

electron. The errors are the independent statistical errors with respect to the nominal fitted mass. 
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